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Abstract 

A brief introduction of the relationship of string percolation to the 
Quantum Chrorno Dynamics (QCD) phase diagram is presented. The behav¬ 
ior of the Polyakov loop close to the critical temperature is studied in terms of 
the color fields inside the clusters of overlapping strings, which are produced 
in high energy hadronic collisions. The non-Abelian nature of the color fields 
implies an enhancement of the transverse momentum and a suppression of 
the multiplicities relative to the non overlapping case. The prediction of this 
framework are compared with experimental results from the SPS, RHIC and 
LHC for pp and AA collisions. Rapidity distributions, probability distribu¬ 
tions of transverse momentum and multiplicities, Bose-Einstein correlations, 
elliptic flow and ridge structures are used to evaluate these comparison. 

The thermodynamical quantities, the temperature, and energy density de¬ 
rived from RHIC and LHC data and Color String Percolation Model (CSPM) 
are used to obtain the shear viscosity to entropy density ratio (rj/s). It 
was observed that the inverse of (77/s) represents the trace anomaly A = 
(e — 3 P)/T 4 . Thus the percolation approach within CSPM can be success¬ 
fully used to describe the initial stages in high energy heavy ion collisions 


* Corresponding author 

Email address: brijesh@purdue.edu (B. K. Srivastava) 


January 8, 2015 



in the soft region in high energy heavy ion collisions. The thermodynamical 
quantities, temperature and the equation of state are in agreement with the 
lattice QCD calculations. Thus the clustering of color sources has a clear 
physical basis although it cannot be deduced directly from QCD. 
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1. Introduction 


1.1. Motivation and historical perspective 

Before the discovery of the quarks there was an interest in the behavior 
of matter at high density and/or high temperature [T]. This interest was 
increased with the formulation of QCD and the possibility of distributing 
high energy or high nucleon density over a large volume to temporary restore 
broken symmetries of the physical vacuum creating abnormal dense states of 
nuclear matter [ 2 ], 

Very early after QCD was born, it was pointed out that the asymptotic 
freedom property of QCD implies the existence of a high density matter 
formed by deconhned quarks and gluons [3] and the exponential increasing 
of the spectrum of Hagedorn was connected to the existence of a differ¬ 
ent phase, in which quarks and gluons are deconhned jl]. The thermalized 
phase of quarks and gluons was called Quark Gluon Plasma (QGP) [5] and 
it was realized that the required high density could be achieved in relativis¬ 
tic heavy-ion collisions HE] and several signatures of QGP were proposed. 
Quarkonium suppression |Hj, anomalous excess of photons and jet quenching 
mm were some of them. At this time, it was pointed out the relevance of 
percolation in the study of the phase structure of hadronic matter mm- 

Experimentally there was a large effort to study in laboratory the decon¬ 
finement and chiral symmetry restoration phase transitions and to explore the 
properties of high density matter starting by the AGS and ISR experiments 
and continued at SPS, RffIC and LffC. The SPS accelerator experiments 
already displayed several signatures that hinted at the onset of QGP forma¬ 
tion [13|. The RffIC data show a striking bulk collective elliptic flow, which 
is generated at relatively early times, since otherwise the spatial anisotropy 
could not convert into a momentum spatial anisotropy. The flow pattern 
was consistent with a very low shear viscosity over entropy density ratio 77 /s, 
indicating strongly interacting matter. On the other hand, jet quenching 
phenomena were clearly observed indicating that this deconhned strongly in¬ 
teracting matter was very opaque The above mentioned 

ratio gave rise to an increasing interest on the AdS/CFT correspondence due 
to the result p/s = 1/47T [T9i] . 

The recent LHC experiments [20, | 2 TJ .22] have extended the study of the 
elliptic how to all the harmonics m Eij conhrming that the high quark 
gluon density matter interacts strongly. The collective behavior and the 
ridge structure observed previously at RHIC in Au-Au and Cu-Cu collisions 
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[25, 26], was also observed in pPb [2711251122] and pp interactions [30] at high 
multiplicity at the LHC. The collective flow of pPb and the ridge structure of 
pPb and pp collisions is a challenge to the usual hydrodynamic descriptions. 
On the other hand, the data on quarkonium seems each time to confirm more 
the validity of the combined picture of a sequential melting of the different 
resonances together with recombination of heavy quarks and antiquarks at 
high energy [31, 32] [33]. Also a detailed study of the jet quenching for identi¬ 
fied particles has been performed ra showing interesting features related to 
the loss of coherence of the gluons emitted in the jet due to the high density 
medium. 

From the theoretical side in addition to the hydrodynamic studies the 
Color Glass Condensate(CGC) picture [35], 36j HZ, HU [39j derived directly 
from QCD, is very appealing and gives a reasonable description of most of 
the experimental observables. At first sight, due to the non-Abelian nature 
of QCD, with gluons carrying color charge, the gluon density xG(x } Q ) rises 
rapidly as a function of the decreasing fractional momentum x, or increas¬ 
ing the resolution Q s . So, the gluon showers generate more gluon showers 
producing an exponential avalanche toward small x. As the transverse size 
of the hadron or the nucleus rises slowly at high energies, the number of 
gluons and their density per unit of area and rapidity increase rapidly as 
x decreases. However, there will be fusion of gluons leading to a limited 
transverse density of gluons at some fixed momentum resolution Q s and that 
is gluon saturation US- The low x gluons are closely packed, the distance 
between them being very small, hence the interaction coupling is weak a s <C 
1. Weak couplings systems are adequate to be studied in QCD. In a given 
collision the multiplicity should be proportional to the number of gluons, 
which at the saturation momentum Q s , is [TP, [T2j. 

dN 1 q q 

7 - 7rr\Qs R 

dy a s {Q s ) 

This dense system, called CGC, has a very high occupation number l/a s (Q s ) 
and corresponds to a highly coherent state of strong color fields. The high x 
gluons can be regarded as the sources of the low x gluons. The independence 
on the cutoff, used to separate the high x gluons from the low x ones, gives 
rise to a kind of evolution equation. In the infinite momentum frame, these 
large momentum gluons travel very fast and therefore their time scales are 
Lorentz dilated. 
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Due to the Lorentz contraction the collision of two nuclei can be seen 
as that of two sheets of colored glass where the color held in each point of 
the sheets is randomly directed. Taking these fields as initial conditions, one 
finds that between the sheets longitudinal color electric and magnetic fields 
are formed. The number of these color flux tubes between the two colliding 
nuclei is Q 2 S R 2 forming Glasrna [13], which has been extensively applied to 
compare with experimental data. 

An alternative approach to the CGC is the percolation of strings m 
351 351 1371 38] subject of this report. The percolation of strings is not di¬ 
rectly obtained from QCD and only is QCD inspired. In string percolation, 
multiparticle production is described in terms of color strings stretched be¬ 
tween the partons of the projectile and target. These strings decay into 
q — q and qq — qq pairs and subsequently hadronize producing the observed 
hadrons. Due to the confinement, the color of these strings is confined to a 
small area S i = 7rrQ, with r 0 =0.2 fm in the transverse space. With increas¬ 
ing energy and/or size and centrality of the colliding objects, the number of 
strings grows and the strings start to overlap forming clusters, very similar 
to discs in 2-dimensional percolation theory [39i|. At a given critical density, 
a macroscopic cluster appears crossing the collision surface, which marks 
the percolation phase transition. Therefore the nature of this transition is 
geometrical. 

In string percolation the basic ingredients are the strings, and it is nec¬ 
essary to know their number, rapidity extension, fragmentation and number 
distribution. All that requires a model and therefore string percolation is 
model dependent. However most of the QCD inspired models give very sim¬ 
ilar results for most of the observables in such a way that the predictions are 
by a large measure independent of the model used. 

The string percolation and the Glasrna are related to each other [ED] . In 
the limit of high density there is a correspondence between physical quantities 
of both approaches. The number of color flux tubes in the Glasrna picture, 
Q;R 2 , has the same dependence on the energy and on the centrality of the 
collision that the number of effective clusters of strings in string percolation. 
In both approaches for the multiplicity distribution a negative binomial dis¬ 
tribution is obtained, where the parameter k that controls the width has the 
same energy and centrality dependence. The role of the occupation number 
1/ay in CGC is played by the fraction of the total available surface covered 
by the strings formed in the collision. The randomness of the color held in 
the CGC gives rise to a reduction of the expected multiplicity. Similarly the 
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randomness of the direction in color space of the color held of the n strings 
of the cluster stipulates that the strength of the resulting held is not n times 
the strength of the individual color held but a /n. This reduction also implies 
a reduction of the multiplicity of particles produced in the collision. Due to 
all these similarities in both approaches the predictions for different physical 
observables are very similar. The string percolation is able to explore also 
the region where the limit of high density has not been reached. 

The plan of the review is as follows: After this introduction in the two 
next chapters the phase structure of QCD and the main results of lattice QCD 
are briehy studied. Percolation studies related to QCD are briehy reviewed 
and later the color string percolation model is introduced. In the follow¬ 
ing sections the applications of the model to the description of transverse 
momentum distributions, rapidity and multiplicity distributions, azimuthal 
distributions, harmonic momenta and ridge structures, Bose-Einstein corre¬ 
lations, J/T suppression, and forward-backward correlations are described. 
Finally, the thermodynamics in the production of particles in string perco¬ 
lation is studied starting by relating the critical percolation density with the 
critical temperature. Afterwards we study the equation of state, the behavior 
of the trace anomaly and ratio between the shear viscosity and the entropy 
density with temperature. 

1.2. Phase diagram of nuclear matter 

What is the behavior of matter when we increase its density? The ob¬ 
served densities of our world expand over many orders of magnitude, from 
~ 10~ 6 nucleons/cm 3 in average in the Universe to ~ 10 38 nucleons/cm 3 
inside a nucleus and ~ 10 39 nucleons/cm 3 in a neutron star. The study of 
the high density limit, specifically the confinement transition from hadrons 
to quarks and gluons can be regarded as the place where the high energy 
collision (two body) probes the short distance limit and therefore meets the 
thermodynamics (many body) of this short distance dynamics [51]. 

Let us start by studying the behavior of a high density medium under 
electromagnetic forces. In a high density medium, the Coulomb potential 
between two electric charges 

V (r) « - (2) 

r 

becomes screened in such a way that the new potential 

e 2 

V(r)^ —exp(— r/r 0 (p)) (3) 

r 


7 


disappears for r > r 0 , where r 0 is the Debye screening radius, which decreases 
as the charge density p increases. This potential change occurs because 
other charges in the medium are around the two original ones, reducing the 
interaction. In this way the interaction range becomes shorter. If we consider 
a bound electromagnetic state, like a hydrogen atom and increase the density 
the Debye radius eventually becomes smaller than the binding radius of the 
atom ro < r and the force between the proton and electron is reduced so that 
the two particles can no longer be bound. An insulator material consisting of 
bound electrical charges at high density becomes a conductor, it undergoes 
a phase transition, the Mott transition (52]. The charge screening dissolves 
the binding of the constituents giving rise to a new state of the matter, the 
plasma. In this transition, the electrons change their mass from the usual 
physical mass in the insulator phase to an effective medium electron mass. In 
fact, the electrons in the conducting phase move in the periodic held of the 
ions, which is very different from the vacuum. In the case of QCD, we have 
confinement, which means that the lines of the color held between two color 
charges are no longer expanded into all the space like in the case of Coulomb 
potential between electrical charges but are confined to a rather small space 
region, forming hadrons, which are colorless bound states of colored quarks. 
However, at high color charge density we expect that color screening takes 
place resulting in an exponential damping of the potential and removing all 
the long range effects of the confining potential. In this way, we expect that 
color screening transforms a color insulator into a color conductor, leading to 
deconhnement of quarks from hadrons 0. This transition, like the insulator 
to conductor in the electromagnetic case, is a collective effect, in which many 
charges participate, and we expect a behavior like a phase transition. 

In a similar way to the electron mass, the quark masses are expected to 
change between the two phases. In the confined phase the quarks are dressed 
with gluons, acquiring an effective constituent mass of around 300 MeV (1/3 
of the proton mass). For light quarks this constituent mass is much larger 
than their mass appearing in the QCD Lagrangian where it is close to zero. 
In another words, the light quark mass in the confined phase is generated by 
the confinement interaction, and when deconhnement occurs this additional 
mass disappears. 

The role of gluons would be to generate the effective quark mass, main¬ 
taining spontaneous chiral symmetry breaking. There will be need of a higher 
temperature T or baryonic potential p to evaporate or melt the gluonic den¬ 
sity of quarks leading to deconkned quarks and gluons with restored chiral 
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Figure 1: Phase diagram of the nuclear matter. Temperature vs baryonic chemical poten¬ 
tial 11551. 


symmetry. According to this picture, from the hadronic phase, first it is 
the dressed quark phase [55] and after the deconfined and the restoration of 
chiral symmetry. In this picture, there are two scales: First, the hadronic 
scale Rh = 1 fm when the density p is below 


1 


Ph 


(4/3)jr Ri 


and the medium is of hadronic nature. For densities p 

1 

Ph< P< Pq 


(4/3)vri?3 


(4) 


(5) 


with R q r^j 


0.3 fm the medium consists of deconfined massive quarks. Finally, 
for densities p > p q the medium becomes quarks and gluons phase with 
restoration of chiral symmetry. 

This picture could be true if the scale R q were independent of T and /!. 
However, the lattice studies have shown, that near /i = 0, color deconfinement 
and chiral symmetry restoration coincide. Hence, in a medium of low baryon 
density, the mass of the constituent quark vanishes at the deconfinement 
point T c and also the screening radius of the gluon cloud size vanishes. At 
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low T and high //, there is no reason to expect a similar behavior and probably 
there will be an intermediate region of massive dressed quarks between the 
hadronic phase and the deconhned and chiral symmetry restoration phase. 
At low T and high /i, other possibilities could exist as quarkyonic matter and 
color superconductivity [5B, EEL SB] • A possible phase diagram is shown in 
Fig. [l] [M] - 

1.3. Lattice QCD 

In finite temperature lattice QCD, the deconfinement order parameter is 
provided by the vacuum expectation value of the Polyakov loop L(x ) defined 
in Euclidean space 

N t 

L(x) = Tr Yl A 4 (x,t) (6) 

t =1 

L(x ) is the ordered product of the SU(3) temporal gauge variables A 4 (x,t) 
at a fixed spatial position, where N t , is the number of lattice points in time 
direction and Tr denotes the trace over color indices. The Polyakov loop 
corresponds to a static quark source and its vacuum expectation value is 
related to the free energy F q of a single quark 



Below the critical temperature T c quarks are confined and F q is infinite im¬ 
plying (L(x))= 0. In a deconhned medium color screening among the gluons 
makes F q finite, hence for T > T c (L(x)) ^ 0. The breaking of chiral sym¬ 
metry is controlled by the chiral condensate X(T) = (TT) ~ M q which 
measures the constituent quark masses obtained from a Lagrangian with 
massless quarks. At high temperature this mass melts, therefore 

x(T)i[ (8) 

T x dehnes another critical temperature. The corresponding derivatives with 
respect to T of (L) and A", the susceptibilities, have been studied in finite 
temperature lattice QCD at vanishing baryon number, showing a sharp peak 
that dehnes respectively T c and T x . The two temperatures, within errors, 
coincide. Also seen is a sharp transition, i.e. a crossover. The quoted value 
of T c is 155 ±9 MeV [591 EDI El]- 
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Energy densities resulting from lattice QCD are shown in Fig. glcft) 
indicating that even for T > 3T C its values are far from the energy density of 
a free gas of quarks and gluons, namely 

£ = ~^[9g++Qq)]T 4 - ( 9 ) 

where g g , g q and g q are the degeneracy numbers of the gluons, quarks and 
antiquarks. This fact indicates that in a rather broad range of temperatures 
the deconhned phase is a strongly interacting medium. This is better seen 
looking at the interaction measure 

A = (e-3 P)/T 4 . (10) 

The energy momentum tensor trace 

T H = + [1 + ( 11 ) 

is (e —3 P) and even for massless quarks ^ 0 as a consequence of the intro¬ 
duction of the scale in the re-normalization process breaking the conformal 
symmetry (trace anomaly). In Fig. [2] (right) the results of lattice QCD are 
shown. A decreases with T very slowly, even less than 1/T 2 . The deconhned 
medium interacts strongly for a broad range of temperatures above T. 


1-4■ Percolation 

A simple example of percolation is the 2-dimensional continuous percola¬ 
tion, which will be used extensively studying string percolation [321 [62]. Let 
us distribute small discs of area 7T7 q randomly on a large surface, allowing 
overlap between them. As the number of discs increases clusters of overlap¬ 
ping discs start to form. If we regard the discs as small drops of water, how 
many drops are needed to form a puddle crossing the considered surface? 
Given N discs, the disc density is p = N/S where S is the surface area. The 
average cluster size increases with p, and at a certain critical value p c the 
cluster spans the whole surface as shown in Fig. [3] Ea . 

The critical density for the onset of continuum percolation has been de¬ 
termined by numerical and Monte-Carlo simulations for different systems, 
which in 2- dimensional case gives 

1.13 
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Figure 2: The energy density and the pressure as function of temperature (left). The 
energy density shows a sharp rise in the temperature region 170-200 MeV. The interaction 
measure calculated using different staggered fermion actions (right) [511 . 


In the thermodynamical limit corresponding to N —> oo, keeping p hxed, the 
distribution of overlaps of the discs is Poissonian with a mean value £ = pnr^, 
£ being a dimensionless quantity 

£ n 

P n = ^rexp(-£). (13) 

n\ 

Hence the fraction of the total area covered by discs is 1 — exp(—£) [47]. For 
the critical value of 1.13 around 2/3 of the area is covered by discs. This 
number 1.13 is obtained for the case of an homogeneous surface [62] 03] • I n 
cases of non homogeneous surface profiles, this factor changes. For example, 
in the cases of circular surfaces with Gaussian or Wood-Saxon profiles the 
critical percolation is reached at 


Pc 


1.5 
7rr o' 


(14) 


Also the fraction of the area covered by strings is no longer given by 1 — 
exp(—£), but by the function [63J 


1 

1 + aexp(—(£ — £ c )/6) 


(15) 
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Isolated Disks 
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Percolation 


Figure 3: Left panel: Disconnected discs, Middle: Cluster formation, Right panel: Over¬ 
lapping discs forming a cluster of communication ED- 


where £ c =1.5. The parameters a and b depend on the profile function. In 
particular b controls the ratio between the width of the border of the profile 
(27T R) and the total area (nR 2 ), therefore is proportional to 1 /R. 

Three-dimensional percolation has been applied to study the phase bound¬ 
aries of high density matter. As mentioned earlier, even before the discovery 
of the quarks, Pomeranchuk [Tj realized that above a certain high density 
hadrons lost their identity. In fact, when the density of a gas of hadrons is 
increased by raising either the temperature or the baryon density, a quark 
of a given hadron will be closer to some quark or antiquark of other hadrons 
than to the original partners. The identity of hadronic matter is lost and now 
there is deconhned quark and antiquark matter. At fi = 0, the percolating 
density for mesons and low density baryons is 

1 24 

* " (iTsFflf" a6/m - (16) 


where Rh= 0.8 fin is the hadron radius. This is the density of the percolating 
cluster at the onset of percolation. If we place overlapping spheres in a large 
volume, the critical density for the percolating spheres is [49| (631; 64j 


Pc 




0.34 

(4/3 )nRi 


~ 0.16 fm 3 . 


(17) 


However at this point only 29% of the space is covered by the overlapping 
spheres and 71% remain empty, very different from that in 2-dimensional 
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percolation. Here, both spheres and empty space form infinite connected 
networks or clusters. The density given by Eq. ( |T7| ) gives the normal nuclear 
matter density. The more accurate critical percolation density for the onset of 
the deconfinement transition is given by Eq. ( fl6| , leading to a density 4 times 
larger than nuclear matter. The existence of two percolation thresholds, 
one for the formation of a spanning cluster of spheres and another for the 
disappearance of a spanning vacuum cluster is a general feature of percolation 
in three or more dimensions. 

Assuming that for the density of 0.6-0.8 firT 3 the cluster is formed by an 
ideal gas of all hadrons and resonances, the temperature of such a gas at this 
density is 170-190 MeV [6S] which is not far from the critical temperature 
obtained in lattice QCD of ~ 155 MeV. The critical density Eq. ( fl7| ) implies 
that the average distance between quark and antiquark at the deconhned 
point is d q ~ 1 j p x J 3 ~ 1.2 fm. 

At high /i and T = 0 one should consider percolation of nucleons having 
an impenetrable spherical hard core, of radius R around R^j 2. Each sphere 
defines a volume V = which is not accessible to the center of any other 
sphere. The spheres can only overlap partially and the distance between their 
centers must be larger than R = 2 R. Then we have again two percolation 
thresholds. Numerical studies show that in the case of the spheres forming 
a spanning cluster there is no variation in the value of the critical density, 
however for the case of the vacuum percolation threshold, now we have [65J 


P = 


(4/3)7rflj 


~ 0.93 fm 


-3 


(18) 


The disappearance of the vacuum cluster for hard spheres requires a higher 
density than needed for permeable spheres. At this point, it is interesting to 
explore the percolation of constituents of mass M q ~ 300 MeV and radius 
R q ~ 0.3 fm. In this case the critical density will be 


P 


1.2 

(4/3)77 


~ 3.5/m 3 . 


(19) 


This value is 4 times larger than the critical value for nucleon percolation 
(0.93 fm ~ 3 ) and 22 times the normal nuclear matter density (0.16 fm~ 3 ). 
According to this, massive deconhned quarks exist between the hadronic 
matter and the deconhned quarks and gluons state. In this intermediate state 
the quarks are deconhned, but the gluons are bound into the constituents 
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quarks. In the lower baryon density limit, the quarks bind into nucleons. In 
the higher baryon density limit a connected medium of deconhned quarks 
and gluons is obtained. Also at fixed baryon density and increasing T the 
deconhned quarks and gluons and the restoration of the chiral symmetry is 
found. 

Let us mention that in the SU(3) lattice gauge theory, spatial clusters 
can be identified as those where the local Polyakov loops (L(x)) have values 
close to some element of the center. The elements of the center group Z 3 , are 
a set of three phases [0, 27 t/3, —27t/3] [ 66 J. Below T c , ((L(x)) = 0) the values 
of (L(x)) grouped around these three phases, show three pronounced peaks 
located at the center phases. Above T c , ( L(x )) 7 ^ 0 the distribution changes. 
One of the peaks grows and the other two shrink. A spontaneous symmetry 
breaking occurs, which leads to a non vanishing ( L(x )). Spatial clusters can 
be defined grouping the sites with a very similar value of (L(x)). The weight 
of the largest cluster increases sharply at T = T c as seen in Fig. |4} Also 
the diameter of the largest cluster, that remains constant below T c , starts to 
rise quickly at T = T c being 40 times larger at T = 1.2 T c indicating that the 
cluster percolates. These results do not change using ensembles with different 
lattice spacing. Therefore in the pure gauge SU(3) theory the deconfinement 
transition is a percolation phase transition (of the second order). 

I 11 the case of full QCD [67], the fermions break the center symmetry ex¬ 
plicitly and act as an external magnetic field in the Ising model, which favors 
the phase 0 of the Polyakov loop, although the other two phases (±27 t/ 3) 
remain populated. As in the pure gauge case at T = T c there is a pronounced 
increase of the dominant phase. The explicit breaking of the center symmetry 
leads to a crossover type of transition. This suggests that the chiral symmetry 
restoration phase transition also can be related to a generalized percolation 
phase transition [ 68 ]. These properties of the Polyakov loop, giving rise to 
domain like structure or clusters of deconfined matter could explain its large 
opacity as well as its near ideal fluid properties, being the origin of the ellip¬ 
tic flow [62] • In experimental collisions at high energy, we expect that color 
strings are formed between the projectile and target partons. These color 
field configurations must have a small transverse size due to confinement. In 
this way, the strings are seen as small discs in the total available surface of 
the collisions. As the number of strings grows with energy and centrality 
degree of the collision, these strings would overlap forming clusters which 
eventually percolate. This 2-dimensional percolation and its phenomenolog¬ 
ical consequences in relation to SPS, RHIC and LHC pp, pA and AA data is 
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Figure 4: Histograms for the distributions of the phase </>(*) of the local loops ( L(x )). The 
left plot shows the distribution below T c while distribution in deconfined phase is shown 
in the right plot [©'. 
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the main subject of this review. In this case the critical percolation density 
is given by Eq. (12) or Eq. (14) in case of realistic profiles. 


1.5. String models 

The phenomenology of string percolation takes its main ingredient, the 
strings, from models, even though most of the predictions are not dependent 
on the details of the models. Majority of the models roughly coincide in basic 
postulates such as the number of strings and its dependence on energy and 
centrality, which is taken from the Glauber-Gribov model. Strings models 
can be divided in models with color exchange between projectile and target 
as the Dual Parton Model (DPM) [70j ED [72], Quark Gluon String Model 
(QGSM) jTH], VENUS [72], EPOS [H], DPM JET [75] and models without 
color exchange where the interaction excites the projectile and target pro¬ 
ducing strings between the partons of both. Examples of this kind of models 
are HIJING [76], PYTHIA [77], AMPT [78], HSP (79], and URQMD [79]. 
We will concentrate in models with color exchange, essentially the DPM, 
which are based on the 1/IV C QCD expansion and are in correspondence via 
unitarity with the Gribov Reggeon calculus. The DPM or QGSM have been 
extensively compared to the experimental data of ISR, SPS and FermiLab 
obtaining an overall agreement ra- 
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In DPM or QGSM, the multiplicity distribution dN/dy of pp collisions is 
described by the formation and fragmentation of 2 k strings m as shown in 

Fig. [5(a). 




Figure 5: (a):Two cut Pomeron diagram (four chain) for proton-proton collisions, 
(b): Single Pomeron exchange and its underlying cylindrical topology. This is a dominant 
contribution to proton-proton elastic scattering at high energies [701 . 


dNPP 1 

V- = + + (20) 

ay a 

where N q ~ qq and N^ qq are the inclusive spectra of hadrons produced in the 
strings stretched between a valence diquark of the projectile (target) and a 
quark of the target (projectile) and N q ~ q are the inclusive spectra of the 
strings stretched between sea quarks and antiquarks. The single particle 
distribution of each string can be obtained by folding the momentum distri¬ 
bution of the partons at the end of the string with the fragmentation function 
of the string 

rl r 1 dN qq ~ q 

N? qq (s,y)= / dx 1 dx 2 pk{xi)pk{x2 )—q- (y-A,s s ), (21) 

Jo Jo dy 

Here y/s~s is the invariant mass of the string, s s = sx±x 2 , where x\ and x 2 
are the light-cone momentum fractions of the constituents at the ends of the 
string. A is the rapidity shift necessary to go from the overall pp center of 
mass (CM) frame to the CM of one string, 

A=bog((l) (22) 

^ x 2 
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The momentum distributions used for the valence quarks, sea quarks or anti¬ 
quarks and valence diquarks are x -1 / 2 , x~ x , and x 3 ^ 2 respectively. In general, 
the distribution of 2 k partons in the proton is 

2k 

Pk{x 1 ,X 2 k , ,X 2 ,X 3 ,....X 2 k-l) = C P k x\ 1/2 X 2 l "" ( 23 ) 

1 


where C k can be found by normalizing p k to unity. 

With these momentum distributions the N qq ~ q and N q ~ qq strings are 
long (due to x 3//2 and x^ 1 / 2 behavior of their extremes) centered at a rapidity 
point shifted with respect to the CM. The N q ~ q strings are short, centered 
at the CM (due to the x -1 of their constituents ends). Concerning the 
fragmentation functions different methods have been used even within the 
same model. In string percolation, mostly the Schwinger mechanism or the 
Lund fragmentation is used. In Eq. (20) a k is the cross section for producing 
2k strings resulting from cutting k Pomerons. 

As the Pomeron has the topology of a cylinder its cutting gives rise to 
two strings (Fig. |5](b)). Using the AGK cutting rules [80], the dependence of 
erfc with the energy is given by 




871 -g exp (A y) 
kz 



k > 0, 


(24) 


where 

2gC exp(Ay) 

R 2 + ay 

and g is the coupling of the Pomeron to the proton, 1+A the intercept of 
the trajectory of the Pomeron, a its slope and C a parameter describing 
the inelastic diffractive states (C =1 means only elastic scattering without 
diffractive states). The total cross section is obtained summing over k 


OO OO _ O 

a tot = exp(Ay) jr jr (^) Z_1 - / p[4o + (-l) 1- *^" 1 ^)]. (25) 

k=0l=k,l>0 Z L - 


The rise of dN/dy is mainly due to short strings, whose number grows 
with energy. On the other hand, outside the region of central rapidity there 
is no contribution of these short strings and the rise with energy is much 
slower. Assuming a Poisson distribution for cutting k Pomerons 

Pk(n) = ^-—j-exp(-kN), (26) 

n\ 
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where N is the mean multiplicity produced when cutting one Pomeron, the 
multiplicity distribution is 

Pk(n) = J2 w k p k(n), w k = —. (27) 

k ° 

Very often < n > P(n) is plotted as a function of n/ < n >. When the result 
is independent of energy one has the well-known Koba-Nielsen-Olesen (KNO) 
scaling. KNO scaling is roughly obeyed up to the highest ISR energy (y/s 
= 63 GeV) but it is clearly violated at SPS, Fermilab, RHIC and LHC. The 
origin of KNO violation can be understood in DPM easily. The contributions 
of multistrings diagrams become increasingly important when s increases, 
and since they contribute mostly to high multiplicities they push upwards the 
high multiplicity tail. The increase with s of the multistrings contributions 
is due both to the increase of the invariant mass of short strings and to the 
s-dependence of the weights. Hence, we expect a larger KNO violation at 
central rapidity region where the short strings contribute than in the whole 
rapidity range or close to the ends of the phase space where the short strings 
do not contribute. 

The width of the KNO multiplicity distributions is related to the fluc¬ 
tuations of the number of strings, which also control the forward-backward 
correlations. These correlations can be described by the approximate linear 
expression 

< rib >= a + brif, (28) 

where rif (&) is the number of particles observed in the forward (backward) 
rapidity window and a and b are given by 

< rib >< n f > 2 — < rifUb >< rif > 

< nj > — < rif > 2 

b= <n f n b >-<n f ><n b > _ 

< nj > — < rif > 2 

Usually the forward and backward rapidity intervals are taken separated 
by a central rapidity window \y\ < y c in such a way that the short range 
correlations are eliminated (y c > 0.5). Consider symmetric forward and 
backward intervals, Syp = Sys = Sy. In any multiple scattering model the 
origin of long range correlations is the fluctuations in the number of multiple 
scatterings 1ST, 82[ [83] . Let N strings decay each into // particles on the 
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average. The slope b can be split into short range correlations (SR) and long 
range correlations (LR) 184 j 

b = bsa + bLR = rtc + TTc' (30) 

where j3 is the ratio between the forward-backward variance and the forward- 
forward variance of the distribution of particles produced in a single string: 


DpB _ < l^fl^b >~ <Hf><^b> 
Dff < > — < /if > 2 


(31) 


and 

c = - , D n —< N 2 > — < N > 2 . (32) 

< iv > J-^FF 

For a large rapidity window between the forward and backward intervals, 
there are no long range correlations in a single string, Dfb = 0 and we have 


b 


1 


i <N>Dff 

1 "T D nI i% 


(33) 


Usually it is assumed that the multiplicity distribution of a single string is 
Poissonian, Dff = Rf and b becomes 


b = 


l + < N > 
Dn/J-F 


(34) 


According to Eq. (34), at low energy there are no fluctuations in the number 


of strings, D N = 0 and b = 0. As the energy increases, D N increases as well 
as b. On the other hand, if we fix the multiplicity we eliminate many of the 
string fluctuations and therefore b is smaller. 

The generalization of DPM to nucleus-nucleus collisions is obtained as 
follows |85| . Consider a collision of a nucleus A with nucleus B in a configu¬ 
ration with Na participating nucleons of A, N B participating nucleons of B 
(assume that Na < Nb) and a total number N c of inelastic collisions. In this 
configuration hadrons are produced in 2N C strings (2 strings for each inelastic 
collision). Of these 2Na are stretched between valence quarks and valence 
diquarks ( q„ — qq^ and qq„ — q^). The remaining N B — Na valence quarks 
and diquarks of B have no valence partner of A and have to form 2 N B — 2Na 
strings with sea quarks and antiquarks of A (qf — qq^ and qf — q^). The 
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remaining 2 N c — 2N B strings are formed between sea quarks and antiquarks 
of A and B (q s — q s ) 


dN AB 

dy 


1 

&AB 


[ a NA,N B ,N c 9( N B 

N a ,N b ,N c 


N A )N A {N qq *- q * (y) + N q v- qq *(y))+ 


(N b - N A )(N q *~ q ?(y) + N q *- qq *(y))+ 

(N c - N B )(N q °~ q ?(y) + N q ~‘ A - q ‘(y)) + sym(N A N B )], (35) 

where ^n^,n b ,n c * s the cross section for Nc inelastic nucleon-nucleon colli¬ 
sions involving N A nucleons of A and N B nucleons of B. This cross section 
have been studied extensively [86] [8?J EBB] together with its different approx¬ 
imations needed for its evaluation. The inclusive spectra, as in the pp case, 
are given by a convolution of momentum distribution and fragmentation 
functions. In the case of A=B we have approximately 

dN AA 

-«< N a > (2 N qq - q *{y) + (2 <k> -2)N qs ~ qs (y))+ 

dy 


(< N c > - < N a >)2 < k > N qs ~ qs (y), 


(36) 


where we have introduced the possibility of having k multiple scattering in 
the individual nucleon-nucleon interactions, which was neglected in the Eq. 


(35). Notice that in the term proportional to the number of collisions it is not 


specified if the collision is soft or hard. In fact there are many inelastic soft 
collisions included in this term. Sometimes it is wrongly assumed that the 
term proportional to the number of collisions contains only hard collisions. 
We observe that in the central rapidity region we have 2 Nk strings which 
for high energy and heavy nuclei is a very large number (more than 1500). 
Due to that we expect interaction among them, and therefore they will not 
fragment in an independent way. We study later such interactions. 

In hadron-nucleus interactions in DPM, QGSM or VENUS the multiplic¬ 
ity distribution P(n) can be approximated by a negative binomial distribu¬ 
tion. 


2. Color strings with fusion and percolation 

2.1. Introduction 

As indicated in Section 1.5, multiparticle production at high energies can 
be described in terms of color strings stretched between the projectile and 
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target HU [SI [SI ESI EDI ED- Hadronization of these strings produces the 
observed hadrons. The basic characteristic feature of the color string model 
with fusion and percolation, which is the subject of this review, is that strings 
are provided with a finite area in the transverse space. In terms of gluon color 
held they can be considered as the color flux tubes stretched between the 
colliding partons, which in the transverse space are restricted to a finite disc 
of a given radius, dictated by the confinement. The mechanism of particle 
creation is then similar to the one in the well-known Schwinger mechanism 
of pair creation in a constant electric held covering all the space, except 
that now the space is finite in the transverse plane. Note that pair creation 
actually splits the space filled with the chromoelectric held into two parts, 
each of them attached to one of the initial and one of the created partons. In 
this way the dynamics of the string evolution consists of successive breaking 
into more strings. 

So creation of particles goes via emission of qq pairs in the held of the 
string. From the start it is relevant to mention one of the important prop¬ 
erties of this mechanism. The transverse dimension of the string is a char¬ 
acteristic independent of the form of the distribution of created partons in 
the transverse momentum, in contrast to what one may think considering 
the string itself as a distribution of partons (say gluons). In the latter case 
the average transverse momentum of emitted particles is obviously inverse 
to the transverse dimension of the source. However it is apparently not 
so for the Schwinger mechanism. With the held constant in the transverse 
plane the spectrum of constituents (photons in the QED case) is just the 5- 
function. However the emitted electrons have non-zero transverse momenta 
whose average is determined by the strength of the held (although the total 
transverse momentum of the pair is indeed zero). Likewise in our color string 
picture emitted partons have average transverse momenta determined by the 
strength of the chromoelectric held and do not depend on the transverse 
dimension of the string. 

At low energies for collision of hadrons and nuclei with relatively small 
atomic numbers the fact that strings have finite dimension has no influence 
on the results. In the transverse plane strings are projected as discs at 
large distances from one another and particle creation does not feel their 
interaction. However with growing energy and/or atomic number of colliding 
particles, the number of strings grows. Once strings have a certain nonzero 
dimension in the transverse space they start to overlap forming clusters, very 
much like discs in the 2-dimensional percolation theory. The geometrical 
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behavior of strings in the transverse plane then follows that of percolating 
discs. In particular at a certain critical string density a macroscopic cluster 
appears (infinite in the thermodynamic limit), which marks the percolation 
phase transition [44J 45,06]. 

The percolation theory governs the geometrical pattern of the string clus¬ 
tering. Its observable implications however require introduction of some dy¬ 
namics to describe string interaction, that is, the behavior of a cluster formed 
by several overlapping strings. 

One can study different possibilities. 

A most naive attitude is to assume that nothing happens as strings over¬ 
lap, in other words, they continue to emit particles independently without 
being affected by their overlapping neighbors. This is a scenario of non¬ 
interacting strings, which closely corresponds to original calculations in the 
color strings approach, oriented at comparatively small energies (and num¬ 
bers of strings). This scenario however contradicts the idea that strings are 
areas of the transverse space filled with color held and thus with energy, since 
in the overlapping areas the energy should have grown. 

In another limiting case one may assume that a cluster of several over¬ 
lapping strings behaves as a single string with an appropriately higher color 
held (a string of higher color, or a “color rope” [92]). This fusion scenario 
was proposed by the authors and later realized as a Monte-Carlo algorithm 
nearly decades ago [93? 5HE5j. It predicts lowering of total multiplicities and 
forward-backward correlations (FBC) and also strange baryon enhancement, 
in a reasonable agreement with the known experimental trends. 

However both discussed scenarios are obviously of a limiting sort. In 
a typical situation strings only partially overlap and there is no reason to 
expect them to fuse into a single stringy object, especially if the overlap 
is small. The transverse space occupied by a cluster of overlapping strings 
splits into a number of areas in which different number of strings overlap, 
including areas where no overlapping takes place. In each such area color 
holds coming from the overlapping strings will add together. As a result the 
total cluster area is split in domains with different color held strength. As a 
hrst approximation, neglecting the interaction at the domain frontiers, one 
may assume that emission of qq pairs in the domains proceeds independently, 
governed by the held strength (“the string tension”) in a given domain. This 
picture implies that clustering of strings actually leads to their proliferation, 
rather than fusion, since each particular overlap may be considered as a 
separate string. Evidently newly formed strings differ not only in their colors 
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but also in their transverse areas. As a simple example consider a cluster of 
two partially overlapping strings (Fig. § m- 



Figure 6: Projections of two overlapping strings onto the transverse plane [47] . 

One distinguishes three different regions: regions 1 and 3 where no over¬ 
lapping takes place and the color field remains the same as in a single string, 
and the overlap region 2 with color fields of both strings summed. In our pic¬ 
ture particle production will proceed independently from these three areas, 
that is, from three different “strings” corresponding to areas 1, 2 and 3. In 
this sense string interaction has split two strings into three of different color, 
area and form in the transverse space. 

We stress that these dynamical assumptions are rather independent of 
the geometrical picture of clusterization. In particular, in each of the scenar¬ 
ios discussed above, at a certain string density there occurs the percolation 
phase transition. However its experimental signatures crucially depend on 
the dynamical contents of string interaction. With no interaction, clustering 
does not change physical observables, so that the geometric percolation will 
not be felt at all. With the interaction between strings turned on, clustering 
(and percolation) lead to well observable implications. 

In this chapter we shall review these implications for the most immediate 
and important observables, multiplicities and transverse momenta spectra of 
produced particles. 

2.2. Multiplicity and transverse momentum for overlapping strings 

As stated in the Introduction, the central dynamical problem is to find 
how the observables change when several strings form a cluster partially 
overlapping. Let us consider a “simple” string stretched between a quark 
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and antiquark with a transverse area Si. It emits partons with the transverse 
momentum distribution 


io(y,p) 


inda 

dyd 2 p 


m 2 ± (p) 

Ce h 


(37) 


where ti is the tension of the simple string, m\ = m 2 + p 2 and m and p are 
the mass and transverse momentum of the emitted parton. In the following 
we mostly consider emitted pions when we take m — 0. In accordance 
with the Schwinger picture of particle emission we assume that tension t± 
is proportional to the field responsible for emission and thus to the color 
charge at the string ends [921 [9 Ij] . For the simple string stretched between 
the quark and antiquark it is proportional to the quark color charge squared 
Q q. According to Eq. (37) the average transverse momentum squared is 


given by tension tp. (pf)^ = t\ and so is proportional to Qo■ We denote 

It is also 


the multiplicity of produced particles per unit rapidity as /i\. 
proportional to the color charge [92j :96j. 

Now consider two simple strings, of areas S\ in the transverse space each, 
partially overlap in the area S^ (region 2 in Fig. [6j, so that = S\—S^ is 
the area in each string not overlapping with the other. A natural assumption 
seems to be that the average color density q of the simple string in the 
transverse plane is a constant q = Qq/S\. For partially overlapping strings 
the color in each of the two non-overlapping areas will then be 


Q i = «s (1) = <?o(5“)/5 I ). 


(38) 


In the overlap area each string will have color 

Q 2 = qS^ =Q 0 (S^/S 1 ). (39) 


The total color in the overlap area Q 2 will be a vector sum of the two over¬ 
lapping colors qS- 2 - In this summation the total color squared should be 
conserved [92]. Thus Q\ = (Q 0 „ + QoJ 2 where Q 0 „ and Q' ov are the two 
vector colors in the overlap area. Since the colors in the two strings may 
generally be oriented in an arbitrary manner respective to one another, the 
average of is zero. Then Q\ = Q 2 0V + Q' ov 2 , which leads to 

Q 2 = y/2qS<® = n/2Q 0 (S ,(2) /5i). (40) 


One observes that, due to its vector nature, the color in the overlap is less 
than the sum of the two overlapping colors. This phenomenon was first 
mentioned in [92] for the so-called color ropes. 
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As mentioned, the simplest observables, the multiplicity /i and the average 
transverse momentum squared (pt')i are directly related to the field strength 
in the string and thus to its generating color. They are both proportional 
to the color [92] [96] . Thus, assuming independent emission from the three 
regions 1, 2, and 3 in Fig. [6] we get for the multiplicity 

/V/A = 2(S (1) /Si) + y/2 (S^/Si), (41) 

where p\ is a multiplicity for a single string. To find (pf'j one has to divide 
the total transverse momentum squared of all observed particles by the total 
multiplicity. In this way for our cluster of two strings we obtain 


/ 2 \ // 2 \ 2 {SV>/S 1 ) + y/2y/2{S'n/S{) 

2{S^/S 1 ) + V2(SW/S 1 ) 


~ 2{SW/S 1 ) + V2 (SM/Si)’ K J 

where {pl^j l is the average transverse momentum squared for a simple string 

and we have used the evident property 2S' (1 ) + 2 = 2S' 1 in the second 
equality. 

Generalizing to any number N of overlapping strings we find the total 
multiplicity as 

p/l^i = (43) 

i 

where the sum goes over all individual overlaps i of n t strings having areas 
S^K Similarly for the (p'i) we obtain 


2\,( 2 \ = EhM^V^i) N 

Vt / / \ Pt ^ E t (S^/S,) E l (SV/S,) • 


In the second equality we again used an evident identity E,; n % — NSi. 
Note that Eqs. (43) and (44) imply a simple relation between the multiplicity 
and transverse momentum 



(45) 


which evidently has a meaning of conservation of the total transverse mo¬ 
mentum produced. 
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Equations (43) and (44) do not look easy to apply. To calculate the sums 


over i one seems to have to identify all individual overlaps of any number of 
strings with their areas. For a large number of strings the latter may have 
very complicated forms and their analysis presents great calculational diffi¬ 
culties. However one immediately recognizes that such individual tracking 
of overlaps is not at all necessary. One can combine all terms with a given 
number of overlapping strings iii = n into a single term, which sums all such 
overlaps into a total area of exactly n overlapping strings S'* 0 *. Then one 


finds instead of Eq. (43) and (44) 


N 


x/oi = £ (sy/so 


n=1 


and 


/(0 1 = 


N 




(46) 


(47) 


In contrast to individual overlap areas S ^ the total ones S'* 0 * can be easily 
calculated (see next subsection). Let the projections of the strings onto the 
transverse space be distributed uniformly in the total interaction area S' with 
a density p. Introduce a dimensionless parameter (“percolation parameter”) 


£ = pS l = N Si/S. 


(48) 


From Eq. (46) we then find that the multiplicity is damped due to over¬ 


lapping by a factor 


m = 


p 


Np i 


£ 


(49) 


where the average is taken over the Poissonian distribution Eq. (13). 

The behavior of F(£) is shown in Fig. [7](a). It smoothly goes down from 
unity at £ = 0 to values around 0.5 at £ = 4 falling as 1 /a/£ for larger £’s. 
According to Eq. (47) the inverse of F shows the rise of the (/Pt')- 

Note that a crude estimate of F(£) can be done from the overall com¬ 
pression of the string area due to overlapping. The fraction of the total area 
occupied by the strings according to Eq. (13) (see also 122) is given by 


£ V n = i - e-f. 


(50) 


n=1 
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Figure 7: (a). Damping of the multiplicity as a function of £. (b): Percolation dispersion 
squared of the multiplicity per string in units as a function of £ [37]. 


The compression is given by Eq. (50) divided by £. According to our picture 


the multiplicity is damped by the square root of the compression factor, so 
that the damping factor is 


F(0 = 




£ 


(51) 


For all the seeming crudeness of this estimate, Eq. (51) is very close to the 
exact result, as shown in Fig. [7](a) by a dashed curve. 


2.3. Multiplicities and their dispersion 

In our picture in the transverse space simple strings are represented by 
discs of radius r 0 and area Si = 7rrQ homogeneously distributed in the total 
area S'. We normalize S' assuming that centers of the discs are inside the 
unit circle of area So = 7r so that S' = 7r(l + a) 2 , where a is the normalized 
disc radius. The disc density is p — N/S and the percolation parameter is 
£ = pS\ = NSi/S. In the thermodynamic limit N —>■ oo, so that at fixed 
£ the radius of the discs goes to zero. For fixed £, a 2 = £/(VN — i/£) 2 , so 
that with our normalization at large N a 1/y/N and Si ~ 1/N. Since the 
discs are distributed homogeneously, the probability that their centers are at 

28 








points Vi, i = 1 ,...N inside the unit circle is independent of r, and is given 

by 

P(r,) = 1/S„ N . (52) 

With the disc centers at points r*. the overlap area of exactly n discs is given 
by the integral 


N 


S t ° t {r l ,...r N ) = J d 2 r n 6>(a — |r — r ifc |) 0 (l r - r *J- a )- 

k=l k=n -\-1 

(53) 

The average of S^ ot will be given by a multiple integral over r* with the 
probability Eq. (52): 


1 N 

( S «‘) = W n dVrtr-i, •■■») = C% / d 2 rF*(r)(l-F(r)) N -\ (54) 

Oq JSo i=1 JS 

where 

F(r) = (l/S'o) / d 2 ri0(a — \r — vA). (55) 

J So 

The function S'o-F(r) gives an area occupied by a circle C of radius a with 
a center at r which is inside the unit circle So. If r < 1 — a then C is always 
inside So so that 

F{r) = cr 0 /S 0 , 0 < r < 1 — a. (56) 

However for r > 1 — a a part of C turns out to be outside the unit circle, and 


F(r) = <7(r)/S 0 , 1 — a<r<l + a, (57) 

where a(r) < Si is the overlap of the two discs C and S 0 . 

Generally, the overlap of two circles of radii ri and r -2 with a distance r 
between their centers is given by 

a(ri,r 2 ,r) = (l/2)r^(ai -sincxi) + (1 /2)rl(a 2 -sina 2 ), (58) 


where 


1 

2r i 


r 2 _ T 

/ 1 / O 


cos(«i/2) = — [r + —- 2 , 


2 _ r 2 ' 

1 '2 


cos(a 2 /2) = - 1 - (r T ' r 
2r 2 V r 


( 59 ) 


The function a(r) in Eq. (57) is just cr(l,a,r). 
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Equations (54) 
any finite value o: 


-(59) allow to calculate numerically the average 
ivwithout much difficulty. 


5'*°^ for 


In the thermodynamic limit, N —» oo with £ being fixed, the calculation 
of ( S n ) becomes trivial. Indeed then one can neglect the part of integration 
in r with r > 1 — a altogether, with an error \ j CL 1/VN. With the same 
precision one then finds 


(ST) = SC n N (a 0 /Sr(l - a 0 /S) 


N-n 


(60) 


where we have put Sq ~ S. The physically relevant values of n remain finite 
as IV —> oo. So we can approximately take 

C n N = N n /n \, (1 - a 0 /S) N ~ n = ex P (-Na 0 /S). (61) 


We then fold that in the thermodynamic limit the distribution of overlaps in 
n is Poissonian with the mean value given by £ (Eq. ©• 

Calculation of the multiplicity dispersion requires knowledge of the aver¬ 
age of its square. With the centers of the discs at r 1} ... r n it has the form 


p 2 (r 1; ...r N ) = (l/u 2 )(^ y/nS^iri,...r N )) 2 , (62) 

n 


where <S'* ot (r 1 , ...rjv) is given by Eq. (53). Taking the average over the discs 
centers positions we now come to a double integral in r and r' 



n N 

y, i r - r *j) n ^(i r " r *j- a ) 

{u,-in}C{l,...iV} k= 1 k=n+ 1 

m N 

Y n%-|r- r ,l) n °(\ r - r ji\- a )- (63) 

{jl i=l l=m- J-l 

This complicated expression, however, continues to be factorized in all 
r t and can be substantially simplihed. Leaving the details to the original 
derivation in m we present here the final expression in the form of the sum 

i^ 2 ) = 4 Y \l(n + p)(n+p) 

^0 n,m,p 
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where 


C% m ’ p I d 2 rd 2 r' 


r,r')x n (r,r')x m {r , ,r)( 


N—n—m—p, 


r, r'). 


( 64 ) 


y^iTl,m,p _ 

^ AT 


N\ 


n\m\p\(N — n — m — p)\ 
x (r,r') = F(r ) - 0(r,r'), 

£(r, r') = 1 — F(r) — F(r') + </>(r, r'), 

0(r,r') = (l/5 0 ) f 9{a — |r — r 1 |)6 ) (a — |r' — r 1 1 ), 
Js 0 


(65) 


with F(r) defined before by Eq. (55). 

This expression is exact and may serve as a basis for the calculation of 
the average square of the multiplicity at finite N. However the new function 
(j) becomes very complicated when both variables r and r' are greater than 


1 — a. For this reason rather than analyze the general expression Eq. (64) for 


finite N we shall immediately take the thermodynamic limit N —> oo. We are 
in fact interested in the dispersion, not in the average square of multiplicity. 
It is important, since the leading terms in N cancel in the dispersion. So we 
shall study the difference 


D 2 = U 2 )-U 


( 66 ) 


in the limit N —y oo, £ finite. As we shall see, although both terms in the 
right-hand side of of Eq. (66) behave as N 2 separately, their difference grows 
only as N. 


Separating from Eq. (64) the term with p — 0 and combining it with the 


second term on the right-hand side of Eq. (66) we present the total dispersion 
squared as a sum of two terms 


D 2 = D\ + D 


2 > 


(67) 


where 


D\ = f d 2 rd 2 r'[C% m x n (r,r')x m (r',r)( N n m (r,r’)- 

^0 n,m "'•5 


C%C%F n (r)F M (r'){l - F(r)) N ~ n { 1 - F(r')) m ] 


( 68 ) 


and D 2 is given by Eq. (64) with a restriction p > 1. 
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Leaving again the details to the original calculations in m we present 
the final results in the thermodynamic limit. 

The first part of the dispersion squared is 


D\/N = 



2 


+ Te 


-2? 


]T J q 2 RdR[( 1 - A (R)) n+m e^ R) - 1], 


f w n ' m 


(69) 


In the first term the averages are to be taken over the Poissonian distribution. 

(70) 

The second one can be easily evaluated numerically. The second part is 


1 R 

A (R) = —(a — sin a), a = 2arccos—. 
7T 2 


Dl/N 


2 e - 2 « v v- y/(» + p)(m + p) r 
i n\m\p\ 


x [ 2 RdRX p {R)(l- X(R)) n+m e^ R) . (71) 

Jo 

This part is evidently positive. Its numerical evaluation shows that it nearly 
cancels the large negative contributions from D\. So the numerical calcula¬ 
tion of the dispersion requires some care. 

As mentioned in Sec. 1.4 percolation is a purely classical mechanism. 
Overlapping strings form clusters. At some critical value of the parameter £ 
a phase transition of the 2nd order occurs: a cluster appears which extends 
over the whole surface (an infinite cluster in the thermodynamic limit). The 
critical value of £ is found to be £ c ~ 1.12 — 1.20 02). Below the phase 
transition point, for £ < £ c , there is no infinite cluster. Above the transition 
point, at £ > £ c an infinite cluster appears with a probability 


l\ = #« - «(? - (of- 


(72) 


The critical exponent (3 can be calculated from Monte-Carlo simulations. 
However the universality of critical behavior, that is, its independence of 
the percolating substrate, allows to borrow its value from lattice percolation, 
where f3 = 5/36. 
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Cluster configuration can be characterized by the occupation numbers 
(z/ n ), that is average numbers of clusters made of n strings. Their behavior 
at all values of £ and n is not known. From scaling considerations in the 
vicinity of the phase transition it has been found EH 


(v„) = n T F(n”((, - y), |£ - y << 1, n » 1, 


(73) 


where r = 187/91 and a = 36/91 and the function F(z ) is finite at z = 0 and 
falls off exponentially for \z\ —> oo. Eq. (73) is of limited value, since near 


£ = £ c the bulk of the contribution is still due to low values of n, for which Eq. 


(50) is not valid. However from Eq. (73) one can End non-analytic parts of 


other quantities of interest at the transition point. In particular, one Ends a 
singular part of the total number of clusters M = as A (M) = c|£—£ c | 8 / 3 . 
This singularity is quite weak: not only (M) itself but also its two first 
derivatives in £ stay continuous at £ = £ c and only the third blows up as 
|£ — £ c | -1 / 3 . So one should not expect that the percolation phase transition 
will be clearly reflected in some peculiar behavior of standard observables. 

Indeed we observed that neither the total multiplicity nor (pf'j show any 
irregularity in the vicinity of the phase transition, that is, at £ around unity. 
This is not surprising since both quantities reflect the overlap structure rather 
than the cluster one. The connectedness property implied in the latter has 
no effect on these global observables. 

It is remarkable, however, that the fluctuations of these observables carry 
some information about the phase transition. As discussed above, the dis¬ 
persion of the multiplicity due to overlapping and clustering can easily be 
calculated in the thermodynamic limit. The result is shown in Fig. [7](b). 
The dispersion shows a clear maximum around £ = 1 (in fact at £ ~ 0.7). 
This value is somewhat lower than the critical one £ ~ 1.2 but still conveys 
certain information about the percolation phase transition around this point 
in spite of the fact that it basically does not feel the connectedness properties 


of the formed clusters. Of course, due to relation Eq. (47), the dispersion of 
■)> has a similar behavior. 

We have to warn against a simplistic interpretation of this result. The 
dispersion shown in Fig. [7](b) is only part of the total one, which besides 
includes contributions from the fluctuations inside the strings and also in 
their number. Below we shall discuss the relevance and magnitude of these 
extra contributions. 

An intriguing question is a relation between the percolation and formation 
of the quark-gluon plasma. Formally these phenomena are different. Perco- 
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lation is related to the connectedness property of the strings. The (cold) 
quark-gluon plasma formation is related to the density of the produced par¬ 
ticles (or, equivalently, the density of their transverse energy). However in 
practice percolation and plasma formation go together. In fact, the trans¬ 
verse energy density inside a single string seems to be sufficient for the plasma 
formation. Percolation makes the total area occupied by strings compara¬ 
ble to the total interaction area, thus, creating, a sizable area with energy 
densities above the plasma formation threshold. 

Let us make some crude estimates. Comparison with the observed mul¬ 
tiplicity densities in pp(p) collisions at present energies fix the number of 
produced (charged) particles per string per unit rapidity at approximately 
unity. Taking the average energy of each particle as 0.4 GeV (which is cer¬ 
tainly a lower bound), formation length in the Bjorken formula [7] as 1 fm 
and the string transverse radius as 0.2 fm [43] we get the transverse energy 
density inside the string as ~ 3 GeV /fin 3 . The plasma threshold is currently 
estimated to be at 1 GeV/ fm 3 . So it is tempting to say that the plasma 
already exists inside strings. This however has little physical sense because 
a very small area is occupied by a string. One can speak of a plasma only 
when the total area occupied by a cluster of strings reaches a sizable fraction 
of the total interaction area. In Fig. [8] we show this fraction for a maximal 
cluster as a function of £ calculated by Monte -Carlo simulations in a system 
of 50 strings. It grows with £ and the fastest growth occurs precisely in the 
region of the percolation phase transition: as £ grows from 0.8 to 1.2 the 
fraction grows from 0.3 to 0.6. With a string cluster occupying more than 
half of the interacting area, one can safely speak of a plasma formed in that 
area. 

2-4- Distribution in the transverse momentum and quenching 

Clusters of strings may take quite complicated forms in the transverse 
plane varying from a disc corresponding to the simple string to long chains of 
such discs. The question arises what the distribution in transverse momenta 
of emitted partons will be from a cluster. As we have seen in Subsection 
2.2 the average transverse momentum from a cluster can be determined in a 
comparatively simple manner, especially in the thermodynamic limit: 

(p 2 t ) = (rfh/no- ( 74 ) 

However this does not fix the distribution in p t uniquely. 
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Figure 8: Fraction of the total interaction area covered by the maximal cluster as a function 

of C 071- 


Let us return to our simple picture of the fusion of two simple strings Fig. 
|6j We assume that partons are emitted independently from the three areas, 
corresponding to the overlap and two remaining areas. Each of these areas 
can be visualized as a set of strings of elementary transverse area d 2 r. The 
total transverse momentum distribution from these three areas will be given 
by the sum of the three integrals 


dp tot 

dyd 2 p 






+ 



d?re 



= C{2S w e~~^ + V2S {2) e * 1 ^} 

and the distribution of emitted particles from the cluster will be 


i{y,p) 


m m j_ 

1 dp‘° l _ 2Sl'>e“Tr + s/2S^e~^ 
ji (lyd 2 p 2SW + V2S(2) 


(75) 


(76) 


Generalizing to many clusters made of different number of simple strings 
we find the transverse momentum distribution in the general case. 



s„.i yH(sy/s 1 )e- m i /(1 -^ i > 

E„=i MsS*/Sy) 


(77) 
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We recall that S * ot is the total area of overlaps of n simple strings. 

We stress that the distribution from the clusters remains isotropic in the 
transverse space in spite of the fact that clusters themselves have different 
forms and their distribution may not be isotropic at all. 

Equation (77) may be used in the Monte-Carlo simulation. For many 


practical problems it can be calculated in the thermodynamic limit 


Pt) = C 


where averaging is done with the Poissonian distribution Eq. 
reasonable accuracy it may be further simplified to 


(78) 


(13). With 


pi 


= Ce 


y/n 


(79) 


which implies that the distribution has the same Gaussian form as for a 
simple string with appropriately enhanced tension. 

Experimental data indicate however that the transverse momentum spec¬ 
tra of emitted particles are not isotropic. Their dependence on azimuthal 
angle together with the anisotropy of the string distribution leads to the 
well-known azimuthal flows (Section 3.3). In view of this fact our string pic¬ 
ture requires certain refinement. As a source of azimuthal anisotropy one 
may introduce quenching of produced partons in the external chromoelectric 
held created by strings. 

Turn again to fusion of two strings (Fig. [6]). Let the observed parton be 
emitted in different azimuthal directions either from the overlap or from the 
remaining part of one of the strings (Fig. [9]). It is clear that the emitted 
partons have to travel paths of different longitude before they go out and 
are observed. Besides, partons going through the overlap meet stronger held 
than those going only through the held of simple strings. So one concludes 
that if partons loose their energy passing through the held their observed 
distribution will depend on their azimuthal angle although initially they were 
emitted isotropically. This implies that the distribution in the transverse 
momentum has the form Eq. (37) with p —> po where Po(p) is the parton 


momentum at the instant of its creation and p is the observed momentum 


P(p) = Ce 


PotpV 


(80) 
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Figure 9: Emission of particles in different azimuthal directions from two overlapping 
strings. 

The dependence of p 0 (p) on the observed momentum will depend both on the 
longitude of the path traveled by the parton and on the strength of the field 
it meets along this path. The concrete form of this dependence is determined 
by the mechanism of quenching. 

Radiative energy loss has been extensively studied for a parton passing 
through the nucleus or quark-gluon plasma as a result of multiple collisions 
with the medium scattering centers [1551 [55] , In our case the situation is 
somewhat different: the created parton moves in the external gluon field 
inside the string. In the crude approximation this field can be taken as being 
constant and orthogonal to the direction of the parton propagation. In the 
same spirit as taken for the mechanism of pair creation, one may assume that 
the reaction force due to radiation is similar to the one in the QED when a 
charged particle is moving in the external electromagnetic field. This force 
causes a loss of energy, which for an ultra-relativistic particle is proportional 
to [its momentum xfield] 2 / 3 [100] : 

^ = - 0.12 eHeEp^f 3 , ( 81 ) 

ax v J 
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where E is the external electric held. Eq. (81) leads to the quenching formula 


Po(p, 0 = p( 1 + i^P 1/3 t 2/3 /) 3 , (82) 

where we identihed eE/w = t as the string tension and l is the longitude of 
the path traveled by the parton in the held with tension t. The quenching 
coefficient k has to be adjusted to the experimental data. In our practical 
calculations it was chosen to give the experimental value for coefficient v 2 
in mid-central Au-Au collisions at 200 GeV, integrated over the transverse 
momenta. 

Of course the possibility to use electrodynamics formulas for the chromo¬ 
dynamic case may raise certain doubts. However in mm it was found that 
at least in the iV = 4 SUSY Yang-Mills case the loss of energy of a colored 
charge moving in the external chromodynamic held was given by essentially 
the same expression as in the QED. 

Note that from the moment of particle creation to the moment of its pas¬ 
sage through other strings a certain time elapses depending on the distance 
and particle velocity. During this time strings decay and the traveling par¬ 
ticle will meet another string partially decayed, with a smaller color Q than 
at the moment of its formation. So one has to consider a non-static string 
distribution with string colors evolving in time and gradually diminishing 
until strings disappear altogether. To study the time evolution of strings we 
again turn to the Schwinger mechanism. For it one has the probability of 
pair creation in unit time and unit volume as mm 

IVt = q— t 2 e~ — i (83) 

47T 

where again t stands for eE/ir in QED. For a realistic string the volume 
V = SL Z where S is the string transverse area and L z is the longitudinal 
dimension of the string. For the single string of color Q we have t = Qt\ 
where t\ is the string tension of the ordinary string with Q = 1. The average 
transverse momentum squared of the emitted quark-antiquark pair < Pi > 
is just t. To estimate L z we assume that the string emits a pair when its 
energy is equal to 2 < p 0 >= \/t, which gives L z = 1/y/t, so that we get the 
average probability in unit time 

T r = t 3/2 S . (84) 
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The string color diminishes by unity with each pair production. So we find 
an equation which describes the time evolution of the string color Q(T) 


dQ(T) 

dT 


with the solution 


Q(T) = 


= -aQ 3/2 (T) 

_ Qo _ 

(1 + i«Tv^) 2 ’ 


(85) 


( 86 ) 


where Qo is the initial color at the moment of the string creation. Coeffi¬ 
cient a = t^ 2 S /(27t) depends on the string transverse area S. In practical 
calculations we use the picture in which the fused string is in fact modeled 
by a set of “ministrings” formed at intersections of simple strings with the 
same area as the simple string, but greater color. This gives a = 0.03 fm~ l . 
The average color of ministrings is of the order 2-3. So it changes only by 
30 -50 % even when the emitted parton travels 5 fm of distance. So the 
time scale of string evolution is estimated to be considerably greater than 
time intervals characteristic for partons traveling inside the string matter. 
However the effect of string decay with time is noticeable and we take it into 
account in our calculations. In fact it practically does not change the results 
but changes the value of the quenching coefficient k, which in any case is 
to be adjusted, as explained above. In this sense our results are practically 
independent of the concrete choice of a in the reasonable interval of values. 

Finally we note that the Schwinger formula Eq. (37) describes well the 
spectra only at very soft po- To extend its validity to higher momenta one 
may use the idea that the string tension fluctuates, which transforms the 
Gaussian distribution into the thermal one |1U31 104j : 


_ pq(p) 

Pip) = Ce . 


(87) 


2.5. Rapidity dependence 

Color strings are stretched between partons into which colliding hadrons 
or nuclei pass before collisions. Let CM energy squared of the collision be 
s = 2 pp ~ 2 p+p- where p and p are the 4-momenta of colliding nucleons 
and we neglect here the nucleon mass m at high energies. A string carries 
a fraction of s depending on the momenta of the partons k with k + = x + p\ 
and k with(fc) = X-p, so that its CM energy squared is s s = sx + x-. Defining 
upper and lower rapidities as 

Y _ Y 

y = — +lnx + , y = -—~ In a;- (88) 
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where Y = ln(s/m 2 ) is the overall rapidity one can say that the string is 
stretched in the final interval of rapidity between y and y. If this interval is 
large then the probability of particle emission at rapidity y p from the string 
will be practically independent of rapidity while y < y p < y. This does not 
imply that the observed particle spectrum will be rapidity independent. In 
fact, as explained in Section 1.5, partons forming the string are distributed 
in the colliding hadrons with probabilities p(y) and p(y), which are obtained 
from Eq. (23) upon integration over spectator variables. The final distribu¬ 
tion in rapidity y p will then be governed by the factor 


r Y rvp 

'PiVp) = d y dyp(y)p(y). (89) 

J y P J y 

and is thus totally determined by the distribution of partons in the colliding 
hadrons. 

Now let us consider string fusion. Let two strings be stretched between 
partons from the projectile with momenta k\ and k 2 and target momenta k\ 
and k 2 with k i+ = Xi + p and = Xj_, i — 1, 2. Conservation of momentum 
dictates that if the two strings completely fuse then the ends of the fused 
string have light-cone momenta (# 1 + + x 2+ )k + and (xi_ + x 2 -)k-. In terms 
of rapidities the new ends will be 

yfused = In (e yi + e V2 ^ , y fuse d = ~ In (e~ yi + e~ V2 ^j 

where y r and y*, i = 1,2 are ends of the fusing strings. This result trivially 
generalizes to fusion of arbitrary number of strings. If end rapidities of the 
strings are y* and y* then if they completely fuse the ends of the fused string 

will be 

Vfused bl ^ ) ( G y ^ , yfused bl ^ ) ) 6 y ^ 

i i 

For illustration, if all the fusing strings are the same then 


yfused = y + In n, yfused = y -Inn 

The situation obviously complicates if strings are not fused completely, as 
in Fig. [6j Then one has to consider the three parts with areas S^ 2 \ and 
as independent strings and for each of them determine ends in rapidity 
separately taking into account the sum of the parton momenta in each of the 
three. 
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3. Model Results and comparison with experiments 


3.1. Multiplicity distributions 

The multiplicity distributions in the DPM or QGSM in pp and AA col¬ 
lisions are given by Eqs. (20 21) and Eqs. (35 36) respectively. However, as 
the energy or the centrality of the collision increases one expects interaction 
among the strings stretched between the projectile and target partons. As 
discussed earlier, due to the randomness of the color field in color space, the 
strength of the resulting color field in a cluster of n strings, is only y/n times 
the strength of the color field of a single string, giving rise to a suppression of 
the multiplicity of particles produced in the decay of the cluster. The same 
reason lies at the origin of the enhancement of the mean transverse momen¬ 
tum with string density. The corresponding equations for both quantities, 
were obtained in the previous section, Eq. (49) and Eq. (74) respectively. At 



Figure 10: Pseudorapidity distribution of charged particles for central Pb-Pb events at a 
beam energy of 3 TeV per nucleon |105l . 

present, most of the models have incorporated the suppression of multiplic¬ 
ities compared with the superposition of independent scatterings result by 
means of different mechanisms, but this was not so twenty years ago, when 
the predictions of all the models for central Pb-Pb collisions at LHC energy 
were a factor two higher than the prediction of the string fusion model (SFM) 


41 



















,1m, which was a previous version of the percolation string model, as it 
is seen in Fig. 10 The SFM prediction is close to the LHC data for Pb-Pb 


collisions at 2.76 TeV. According to Eq. (49) the multiplicity distribution in 


pp collisions in the central rapidity region is given by 


dNPP 

dy 


= F (Q N pVu 


( = N s — 

s V p S ' 


(90) 


where factor F is given by Eq. (51), Np is the number of strings in the 
central rapidity region, Si = 7rrg with ro = 0.2 fm and S is the transverse 
area of the proton. In the DPM in nucleus-nucleus collisions the number of 
strings stretched between the sea quark and antiquarks in the central rapidity 
region is proportional to N A 3 — N A and is given by Eq. (36). Here N A 3 is 
the total number of nucleon-nucleon collisions. Hence, the total number of 
strings in a central heavy ion collisions can be very large, actually more than 
~ 10 3 strings are produced in Au-Au at RHIC energies. However, each string 
must have a minimum of energy to be produced and decay subsequently into 
particles (at least two pions). On the other hand, the total energy available 
in a collision grows as A, whereas the number of strings A 4 / 3 in central 
collisions. Hence, at not very high energy (for instance RHIC energy) the 
energy is not sufficient to produce such a large number of strings. In order to 
take into account this energy momentum conservation effect one may reduce 
the number of sea quark and antiquark strings, changing [TUB] 


Na /3 ->• N\ +a{Vs) 


with 



1 


1 + ln(^/(s/so) + 1) 


One can thus write 


dN AA 

dy 




N a (N^- 1 ) 


dNPP 

dy 


(91) 

(92) 


(93) 


Parameter So marks the energy squared below which energy-momentum con¬ 
servation effects become small and a —>1/3. Up to here we do not take into 
account the interaction among the strings. If we do take into account the 
interaction of strings then we can write a closed formula for the multiplicity 
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distribution in AA collisions in terms of the multiplicity distribution in pp 
collisions, namely [106] 


1 dN 


N a dy 


dN pp 


y =o 


dy 


y =o 


( 1 + 


F(fa,) 

HQ 


N' 


a(y's)—1> 


(94) 


where 

6v„ = f p iVy Vi)+1 ^ L (95) 

on a 

and Sjv a is the transverse area of the collision region formed when there are 
Na wounded nucleons of the projectile and Na nucleons of the target. S^ A 
depends on Na and A. The dependence of the multiplicity on the center of 
mass collision energy A fs is fully specified once the average number of strings 
in a pp collision N* is known. At low energy N* is approximately equal to 
2 , growing with energy as (\Zs/m p ) 2X so that 


K = 2 + 4(^-Y(A) 


2A 


Rn 


m r 


(96) 


Here a single parameter A describes the rise of the multiplicity with energy 
for both pp and AA multiplicity distributions, even though in AA central 
collisions the multiplicity increases faster than in pp collisions due to the 
energy dependent factor ex, arising from energy momentum conservation. 

A fit to pp collisions data in the range 53 < a/s < 7000 GeV and to 
AA collisions (Au-Au, Cu-Cu, and Pb-Pb) at different centralities for 19.6 
< \fs < 2760 GeV has been done. The values obtained from the fit for the 
two parameters are a/s 0 = 245 GeV and A =0.201. 

Figure [TT] shows a comparison of the results for the dependence of midra¬ 
pidity multiplicity on the energy with data for pp |1071 108 . Tf)9 . 11T01 Hill 
111211113] and central Cu-Cu {Na = 50, A = 63) |114j and for Au-Au/Pb- 
Pb [Na — 175, A = 200) |115j . In Fig. 12 the result of the dependence 


of the multiplicity per participant nucleon on the number of participants is 
shown together with the experimental data for Cu-Cu (a /snn = 22.4, 62.4, 
200 GeV) for Au-Au (a /s~nn= 19.6, 62.4, 130 and 200 GeV) and for Pb-Pb 
(a /snn = 2.76, 3.2, 3.9, 5.5 TeV). The evolution outside the central rapid¬ 
ity region has been studied extensively extending to all rapidities Eq. (96) 

msumunsung: 


1 dN NANA 
N a drj 


= KJF(QN 


1 + -') 


exp 




+ 1 


(97) 
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Figure 11: Comparison of the evolution of the mid-rapidity multiplicity with energy from 
the CSPM and data for pp and A-A collisions. Lines are from the model for pp (gray), 
Cu-Cu (blue) and red lines for Au-Au/Pb-Pb ll'H'i . 



file 


Figure 12: Multiplicity dependence on centrality ( N part ). Cu-Cu (triangles), Au-Au (stars) 
and Pb-Pb (circles).Curves represent the model calculations. Blue line for Cu-Cu, green 
line for Au-Au and red for Pb-Pb 1106 . 
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where J is the Jacobean 


and 


K = 


J = 

k 

J (r i = 0) 


cosh ri 


\Jk\ + sinli 2 j] 


( exp( 


-(1 -a)Y 


) + l). 


(98) 


Parameters a and <5 are obtained from fitting to the data (a = 0.34, 5 = 0.84). 
The pseudorapidity dependence is described by the same factor 

(99) 


exp (!tM) + 1 

in pp and AA collisions. This dependence gives rise to a smaller increase at 
central pseudorapidity rj — 0 than at large pseudorapidity 7] = Y [ 117. 118j . 
The limiting fragmentation property have been studied carefully, showing it 
is not exact and a violation of it should be more visible at the highest LHC 
energies [ 116 . 117. 118j . In Fig. 13 we show the comparison of Eq. ( |97| ) with 
experimental data at different energies (53 GeV < i/s < 7 TeV) [120. 121j. 


In Fig. 14[a-c) we show the comparison between the results for Cu-Cu [121J. 
Au-Au |122] and Pb-Pb [123] collisions and the experimental data. 

The percolation of strings have been applied to the case of different pro¬ 
jectile and target as well. In Fig. 14 [d) the results for d+Au at different 
centralities compared to the experimental results mi are shown. Similar 
results were obtained in Ref. [124], 

The behavior obtained for dN/dy in pp and AA collisions with energy 
and number of participants is very similar to the Glasma picture in CGC. In 
fact, in CGC the multiplicity distribution is given by 


dN 

dy 


1 

Cls (Q s 


QlH 


■A- 


( 100 ) 


Since the saturation momentum squared Q 2 S behaves like N j 1 / 3 , the multi¬ 
plicity per participant is almost independent of Na and a weak dependence 
arises from the logarithmic dependence on Na of the running coupling con¬ 
stant 1/ols ~ log Na- In percolation the multiplicity per participant is almost 
independent of Na as well and the only additional dependence arises from 
the factor (1 — exp(—^')) 1 A which grows weakly with Na above the perco¬ 
lation threshold. Concerning the energy dependence, Q 2 behaves like s A , so 
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Figure 13: Comparison of results from the evolution of dn c h/dr] with dependence on 
pseudorapidity for pp collisions at different energies (lines) [119] , 


that the same behavior is obtained in percolation. There is an extra energy 
dependence dne to the running coupling constant cc s (Q s ) in the CGC, which 
again corresponds to the energy dependence of the factor (1 — exp(—£)) 1//2 . 
It is not surprising that their exist a correspondence between 1 /a s , the oc¬ 
cupation number or the number of gluons, and the factor 1 — exp(—£) which 
represents the fraction of the collision area covered by strings. The larger the 
occupation number the larger the fraction is. The similarities between the 
Glasma picture of CGC and the percolation of strings are visible not only in 
multiplicities but in most of the other observables, as discussed later. 

In order to explain the faster rise of the multiplicity in central AA colli¬ 
sions than in pp collisions several possibilities have been proposed in CGC, 
such as enhanced parton showers in AA collisions due to the larger aver¬ 
age transverse momentum of the initially produced minijets compared to pp 
collisions |125j or non trivial Q effects intertwined with impact parameter 
dependence [126]. 

3.2. Transverse momentum distributions 

In Section 2.4 we studied the effects of percolation of strings on the mean 
transverse momentum and the dispersion of the transverse momentum dis- 
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tribution and the experimental data related to it. 

As explained in Sec. 2 the detailed study of both multiplicity and mo¬ 
mentum distribution requires analysis of all overlaps of created strings which 
differ both in the number of the strings in the overlap and in the area of a 
particular overlap. With a large number of strings created in AA collisions 
at high energy this is hardly feasible. So it is reasonable to search for some 
effective way to describe the observables in this situation. 






Figure 14: Comparison of results from the evolution of dn c h/dr\ with dependence in pseu¬ 
dorapidity for (a).Cu-Cu at 200 GeV, (b) Au-Au 200 GeV. Plot (c) show ( N 1 1 / 2 ) for 
Pb-Pb collisions at 2.76 TeV (d) dn/d 77 at different centralities for d+Au collisions at 200 
GeV HT51 . 

Let us start with multiplicities. They come from a set of overlaps (“min¬ 
istrings” ) and depend on both the number of overlapped strings and the area 
of the overlap, which combine to give an average multiplicity N from this 
overlap. We may characterize different overlaps just by this average mul¬ 
tiplicity. We call N a “size” of the overlap, the quantity combining both 
the number of overlapped strings and the area. With a lot of overlaps N 
will be changing practically continuously. We then can introduce a proba- 
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bility W(N) to have overlaps with size N in a collision and write the total 
distribution in multiplicity as 

P{n) = J dNW(N)P(N, n), (101) 

where P(N,n ) is the multiplicity distribution from the overlap of a given N, 
which we take to be Poissonian with the average multiplicity N 

P(N,n) = exp(—N)N n /n\. (102) 

The normalization conditions J2 n P( n ) = J2nP(N,n) = 1 and the condition 
Yjn nP(N, n ) = N lead to relations 


J dNW(N) = 1, < n >=< N >- j dNNW(N). 


(103) 


For the weight function we assume the gamma distribution [ 127 , 128, 129j 


W(N) = G(N, k N ,r n ) = exp(-r N N), rjv = 


k 


N 


r (k 


N 


<N> 


(104) 

There are several reason for this choice. First, the gamma distribution re¬ 
produces to a good approximation the cluster size distributions at different 
centralities. In fact, let us consider a peripheral collision, where the density 
of strings is small and there are only very few overlapping strings. In this 
case, the cluster size is peaked at low values of the number of strings of the 
cluster. As the centrality increases, the density of strings increases as well 
,and there are more and more overlapping strings. The cluster size distri¬ 
bution becomes strongly modified. Figure [TS] shows a plot of three different 


cluster size distributions corresponding to three centralities. Each one can 
be described by a gamma function corresponding to different values of kx. 

There is another reason for choosing the gamma function, related to the 
re-normalization group. The growth of the centrality can be seen as a trans¬ 
formation of the cluster size distribution. Start with a set of single strings 
with a few clusters formed of a few overlapping strings. As the centrality 
increases, there appear more strings and more clusters composed of more 
strings. This change can be considered as substitution of strings in a cluster 
by newly formed clusters, defined by new < n >, corresponding to a higher 
color held in the cluster. This transformation, similar to the block transfor¬ 
mations of Wilson type, can be seen as a transformation of the cluster size 
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probability of the type: 


P(x) —> 


xP(x) 

< X >' 


->■ 


x k P(x ) 
< x k > 


->■ ... 


(105) 


Transformations of this kind were studied long time ago by Jona-Lasinio in 
connection with the re-normalization group in probability theory [ 130] , show¬ 
ing that the only probabilities P{x) stable under such transformations are 
the generalized gamma functions. Among the generalized gamma functions 
the simplest one is the gamma function which in addition has one param¬ 
eter less (This additional parameter could be used to refine the model for 
comparison with the data). The transformations of type Eq. (105) have 


been used previously to study the probability associated with some special 
events which are shadowed only by themselves and not for the total of events 
jI5nH3!flH331H331H35] . 

Notice that W(N) satisfies KNO scaling, namely the product NW(N) 
is only a function of N/ < N > and not of the energy. This property is 
a consequence of the invariance of the form of the gamma functions under 


transformations of type Eq. (105) [ 136] . 

Now we pass to the transverse momentum distribution (TMD) f(pt). 
As in the case of multiplicities it conies from overlaps of different number 
of strings having different areas and the TMD from an overlap depends on 
its both characteristics (we may again call this combination a “size” of the 
overlap, this time in respect to the TMD). Take the TMD from an overlap 
of a given size given by the Schwinger mechanism 


fipt,x) = exp (~xp 2 t ) 


(106) 


“size”. 


where x is just this 
can write the total TMD, in all similarity to Eq. (101), as 


Assuming again that x varies continuously one 


f(pt) = / dxW p (x)f(x,p t ) 


(107) 


with a certain positive weight W p (x). To understand its property we use the 
normalization condition for TMD 


This gives a relation 


J dptfipt) =< n > 


dx . . 

— WJx) =< n > . 
x 


(108) 

(109) 
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On the other hand, introducing integration variable N = ax we have from 

( 110 ) 


Eq. <1031 


< n >= j dxxa 2 W(ax). 

Comparing this with Eq. (109[) we can make an identification 


W p (x ) = ( ax) 2 W(ax ) 


( 111 ) 


If we take the gamma distribution Eq. (104) for W(N ) then W p (x) turns out 
up to a factor to be also the gamma distribution bit with different k and r 


W p (x) = —G(x, k + 2, r p ) 


( 112 ) 


with r p = ar (note that a is dimensionful). 

So in the end both the distribution P(ri) and TMD f(p f) are given by a 
convolution of the cluster multiplicity and its TMD with the size probability 
W, which in both cases can be taken as the gamma distribution although with 
different parameter k^,rN and k p ,r p respectively. In the following having in 
mind that most applications will be devoted to the TMD, we denote k p and 
r p as simply k and r leaving notations k! and r' for k n and r n. In many 
discussions the behavior of k and k' is similar and we do not specify which 
of the two fc’s we are discussing. 

Introducing Eq. (104) into Eq. (107) and Eq. (101), we obtain 

| pOO rp 

/ c7rexp(— p 2 x)( )(rx) k ~ 1 exp(— rx) (113) 

Jo T k) 


(1 + ^) 


p ?^k 


and 


T(n + k') 


T(n + l)T(k') (1 + r')(n + k') 


dN 


e N N n 


(r IV)^- 1 ) exp(-r N). 


n\ 


r(fc') 


(114) 


The mean value and the dispersion of the distributions Eq. (113) and Eq 

k < x 2 > — < x > 2 


(114) are 


< x >= 


< n >=< N > —, 
r 


< x > z 
k’ < N 2 > - < N > 5 


1 

k 


< N > 2 


1 

P 


(115) 

(116) 
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( 117 ) 


< n z > - < n > z > 


1 

k' 


<N> 


< n 

The distribution given by Eq. (113) is the negative binomial distribution. 


Equations (113) and (114) can be seen as a superposition of sources (clusters) 
where 1 /k and 1/k' fix the transverse momentum fluctuations and the fluctu¬ 
ations on the number of strings in the clusters. At small string density there 
are almost no strings overlapped, the strings are isolated and k and k' tend to 
infinity. When the density increases, there will be some overlapping strings 
forming clusters, and therefore fc’s decrease. Their minimum is reached when 
fluctuations in the number of strings per cluster reach their maximum. Above 
this point, with increasing string density, these fluctuations decrease and fc’s 
increase. We recall that k and k' are generally different. When comparing 
to the experimental data to fix them, one also has to take into account the 
experimental conditions, such as acceptance and range of variables studied. 

Now we have to take into account that the mean transverse momentum 
and mean multiplicity of a cluster of strings is given by Eqs. (J49l) and (74) 


and therefore we should incorporate this into Eqs. (113) and (114). In this 
way the respective distributions become 


/(*,,) -i 


.2 \ k 


and 


P(n) = 


dpfdy dy k <pf>\ A , ZilALl 

T(n + k') (<n>),F(£)) 


(H8) 


T(n + l)r(fc') ( 1 + /_ w _A n+fe ' 

^ bn>! F{i)>) 


(119) 


We observe that 

< n >= F(£)N S <n>i, < p 2 t >= 


k < pj >j 

k~ 2 W) 


( 120 ) 


Equations (118) and (119) give the distributions for any projectile, target, 


energy and degree of centrality and are universal functions which depend 
only on two parameters <p\>\ and < n>\ the transverse momentum and 
multiplicity of particles produced by one string. In the case of identified 
secondary particles the corresponding quantities for each identified particle, 
< pi >i i and < n >u should be used. In addition to these parameters, 
parameters k and k ', actually are functions of centrality and, as we said 
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Figure 15: Schematic representation of the number of clusters as a function of the number 
of strings of each cluster at three different centralities [128] . 


above, with the growth of string density first decrease up to a minimum close 
to the critical percolation value and afterwards increase . At £ —* oo, k —* oo 
and the distribution given by Eq. (118) becomes exp (—F(£)p 2 / < p 2 > x ) 


very similar to the behavior at £ —» 0. From Eq. (118) we have. 


dlnf 

dlnpt 


-2 F (0 


Pt 


(i , mri ) 

V fc<Pt>ii/ 


< Pt >u 


( 121 ) 


At p 2 —» 0 this reduces to —2 F(^)p 2 / < p 2 >i* and vanishes at p t = 0. On the 
other hand, as the mean p t of particles produced in a single string is larger 
for protons than for kaons and the latter larger than for pions, < p 2 > i7r < 


< Pt >i k < < Pt >i P , the absolute value of the left-hand side of Eq. (121) is 
larger for pions than for kaons and than for protons. This is the well known 
hierarchy that is very often advocated as evidence in favor of hydrodynamics 
models. In Fig. 16 the results of our distribution Eq. (118) are shown 


compared to the RHIC data of PHOBOS [ 137] for central Au-Au collisions 
at yj snn = 200 GeV. We use the values < pf >i n — 0.06, < p 2 >ik— 0.14, 

< p 2 >i p — 0.30 and k = 4. The value of the string density is 2.5. We 
observe a good agreement taking into account that our result is not a fit but 
the result of Eq. (118) fixing the normalization of the three curves at the 
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same value at p t 


2 GeV/c. Let us now discuss the interplay between low 



Figure 16: Experimental PHOBOS data on low p t distributions for pions, kaons and 
protons along with our results for central Au-Au collisions at q snn = 200 GeV [128] . 

and high p t . One defines the ratio Rcp{pt ) between central and peripheral 
collisions as 

AcpM = / f ' (P, ' y A° ) ^°" . (122) 

f(pt,y = o )/N co u 

where the distributions in the numerator and denominator correspond to 
central and peripheral collisions respectively, £ > £. The normalization on 


the number of collisions in Eq. (122), essentially eliminates N s , the number 


of strings, from dN/dy (this is true at midrapidity and not at forward or 
backward rapidities). From Eqs. (|107) and (118) we obtain 


Rcp(Pt) — 


{(k-l)/k) VF(£)J 


A , V 

V 1 k<p$>u) 

O , F^'yp? 

I 1 + k' <p 2 t >u 


(123) 


Here k and k! are values of the parameter k for transverse momentum dis¬ 


tribution for peripheral and central collisions. 
F(0 < F(£) we obtain 


In the limit py —y 0, as 


R cr(0) - (y§y)' 


< 1, (124) 
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which is approximately independent of k and k'. As £'/£ increases the ratio 
Rcp decreases, in agreement with the experimental data. As pt increases, we 
have 

1+ mpi 

~ ylf < 125 > 


k'<pf> li 

and increases with p t (again F(£) > F(£)). At large p t 

d ( \ F (0 2 (k-k') 

RofM ~ ?(F)¥ P - ' 


(126) 


At low density in the region where k decreases with the string density k' < k 
and R-cpipt.) > 1. It is the Cronin effect. As £ /£ increases the ratio Rcp 

111 k + /k + (On¬ 

es 


ar 

k' 




n, 




Figure 17: Ratios for different distributions K/p, p/p in Au-Au collisions at ^/s^n = 200 
GeV at two different centralities :(0-5%, solid circles, 60-70-% open circles) in comparison 
with the data [l29j . 


the normalization of N co u in Eq. (122) does not cancel N s from (IN/dy, 
since in this region N s is proportional to Na instead of N co u- Therefore, an 

additional factor appears now in Rcp{pt)- As N' coll — N' a for central 
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collision is larger than the corresponding quantity for peripheral collisions, 
we have Rcp(pt,y = 3) < Rcp(pt,y = 0), thus a further suppression occurs 
in agreement with the experimental data [138J. 

The results for the TMD for 7 r + , K + and p in Au-Au collisions at — 

200 GeV are in good agreement with the PHENIX experimental data ra. 


In Fig. 17 we show the ratios kaons/pions and proton/pions as a function of 
pt at the two extreme centralities compared to the experimental data. The 
obtained values of k as a function of the string density evaluated for each 
centrality are shown in Fig. 18 The expected increase of k with the string 
density is clearly seen. 

The experimental data on pp in the range yfs =23, 200, 630, 1830, and 
70000 GeV can be described by the distribution Eq. ( |118[ ) m as well. In 


this case the values of k decrease with energy as expected, a higher energy 
means a higher string density. Notice that for string densities above the 
critical percolation value, k should increase, hence in pp a change in the 
behavior of k is nredicted. At 14 TeV the string densitv would be close to 
such critical i 

k«) 



Figure 18: fc(0 extracted from 7r + fits. The behavior expected in percolation is suggested 
by the data (12 9] . 

Concerning the baryon p t distributions, the difference with the meson 
spectrum is not only due to the difference of masses which results in that 
< Pt > Mesons < < Pt >Baryons■ This effect only causes a shift of the maximum 
of the nuclear modified factor Raa but keeps the same height at the maximum 
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contrary to the data. In the fragmentation of a cluster formed by several 
strings the flavor properties follow from the corresponding properties of the 
individual strings and hence the resulting flavor of the cluster F is the flavor 
composition of the individual strings / as well as the color composition. 
We obtain clusters with higher color and different flavored ends. In the 
fragmentation of a cluster we consider the creation of a pair of complexes 
FF. After the decay the two new FF strings are treated in the same way 
and therefore decay into more FF strings until they come to objects with 
masses comparable to hadron masses, which are identified with observable 
hadrons by combining the produced flavor with statistical weights. In this 
way the production of baryons and antibaryons is enhanced. We observe that 
additional quarks or antiquarks required to form a baryon or antibaryon are 
provided by the quarks or antiquarks of the overlapping strings that form 
the cluster. In some sense this coalescence is incorporated in a natural way 
[ 141 1 . 1421 1143] in the CSPM. In order to take this into account, keeping a 
similar formula for the TMD, we modify the latter in an effective way [ 144] . 
We take for the multiplicity (per unit rapidity) of baryons and antibaryons. 


d F = N 

dy 


1 +a 


F(0 < n >i B 


(127) 


instead of Eq. ( |107 ) fitting the parameter a to reproduce the experimental 
data (a =0.09). From Eq. (127) we observe that when a (anti)baryon is 
triggered the effective number of strings is N} +a instead of N s . This means 
that the string density £ must be replaced by 1V"£. The (anti)baryons probe 
a higher density than the mesons for the same energy. The value of a is dif¬ 
ferent for baryons and antibaryons because the rise with centrality is slightly 
different for them, depending on the specific kind of baryon or antibaryon 
as well. In Fig. 19 we plot the results [ 144] for R C p (0-10% central)/(60- 


92% peripheral) for pions and (p + p)/2 compared to PHENIX data. At 5.5 
TeV the LHC predictions are shown. In Fig. 19 we plot the antiproton to 


neutral pion ratio for central and peripheral centrality bins compared to the 
PHENIX data together with the predictions at 5.5 TeV. More results can be 
seen in Ref. [145]. where the the open charm production has been studied. 
Related to percolation, there are other extensive studies of the effects of the 
strong color held on the string tension, including the effect of baryon junc¬ 
tion. These studies extend to strangeness production [1461 1147] and charm 
production [ 148 . 149 ] in pp and AA collisions. 

Equation (118) for the pt distribution does not include the high p t region 
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Figure 19: (a )Rcp (0-10% central/60-92% peripheral) for pions(solid line) and (p+p)/2 
(dashed) compared to the PHENIX data. Blue line is LHC predictions, (b) p/7r° as a 
function of p t for 0-10% (solid) and 60-92% (dashed) centrality bins. The points are from 
PHENIX data. LHC predictions are in blue |144| . 


as far as we take a Gaussian shape for the fragmentation of a cluster. A more 
detailed form for the cluster taking into account perturbative QCD for the 
high p t region may extend the range of validity of our modified distribution. 

There are several scaling properties found in the TMD related to string 
percolation. The experimental data for pp collisions exhibit a universal be¬ 
havior in a suitable variable z = pt/K 


cf)(z) = A 


d 2 N 


2 p t dp t dp 


I pt=Kz- 


(128) 


ffere the parameters A and K depend on the energy of the pp or pp collisions. 

fs =0.9, 


This universal scaling behavior [EDI. as is shown in Fig. f20kor pp at 
2.36 and 7 TeV and pp at 0.63, 1.8, 1.96 TeV. The value of K increases with 
the energy. In string percolation this scaling is satisfied by the distribution 
Eq. (118) with K proportional to 1/F(£). As F(£) decreases with energy K 
increases. A similar scaling has been found in Au-Au (1511 1152) at different 
centralities for the distributions of kinetic energy of pions, kaons and protons. 
Again the percolation of strings distribution satisfies this scaling. 
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Figure 20: Scaling behavior of the charged hadron p t spectra presented in z (a) in pp 
collisions and (b) pp collisions with different energy scales. The inset is the distribution 
of the ratio between the experimental data and the fitted results |150j . 


Recently, it has been shown that in pp collisions the p t spectra of charged 
particles exhibit geometrical scaling, which is also satisfied in nucleus-nucleus 
collisions [1531 1154] . Using proper normalizations it was shown that in any 
type of collision for any centrality and energy the different distributions ap¬ 
proximately are lying on the same curve for low and moderate pt [155] , namely 


1 dN 
N a dpi 





(129) 


This is illustrated in Fig. 21, where the data for pp, pPb, Cu-Cu, Au-Au, 


and Pb-Pb at different energies and centralities are plotted. The percolation 
of strings, at low p t satisfies this scaling. In fact at low p t , the distribution 
Eq. (118) becomes exp(— p t F{^)/ < p t >) which is a function only of Pt/Q s 
if Q s is proportional to < p t > /F 1//2 (£). 

The experimental data on the mean p t as a function of the multiplicity 
have shown an interesting systematic in pp, pPb, and Pb-Pb collisions at 
LHC. All of them grow with multiplicity, being larger in pp than in pPb and 
in Pb-Pb collisions [1561 1157] . On the other hand, in Pb-Pb collisions the 
rise of mean p t with multiplicity is flattened above a certain low multiplicity. 
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Figure 21: Charged particle multiplicity per participant at pseudorapidity 0.2 < 77 < 1.4 
for Au-Au and Cu-Cu central collisions at two RHIC energies 62.4 and 200 GeV. and for 
Pb-Pb collisions at 2.76 TeV plotted as a function of r for A = 0.27 [1551 . 


The same occurs for pPb collisions, even though in this case the change of 
behavior is at higher multiplicity. All this behavior is understood easily in 
string percolation as a consequence of Eq. (120). In fact factor 1/F 2 (£) is 
responsible for the rise with multiplicity of mean p t because £ grows with 
multiplicity. The flattening of the mean p t exhibited in Pb-Pb and pPb data 
is due to the dependence of the parameter k in Eq. (120) on the string density 
£. In Pb-Pb for most of multiplicities the corresponding string densities £ 
are above the critical percolation value. In this region, k grows with £ and 
the rise of p t is lowered. The same happens in pPb, but now we need higher 
multiplicities to have string densities above the critical percolation density. 
In the case of pp collisions most of the corresponding string densities lie 
below the critical percolation string density. In this region k is a decreasing 
function of £, hence there is no flattening. We expect that at higher energies 
even in pp collisions the critical density will be reached. However at 14 TeV 
probably the string density will be slightly below. This discussion refers to 
charged particles but can be extended to identified particles. In the case 
of baryons, as noted before, there are additional enhancements due to the 
higher flavor of the cluster. This can be taken into account in an effective 
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way assuming that baryons probe a higher string density. Due to that we 
expect a flattening of the mean p t of baryons as a function of the multiplicity 
at future LHC energies. 


3.3. Collective flow and ridge 

The cluster formed by the strings has generally an asymmetric form in the 
transverse plane and acquires dimensions comparable to the nuclear overlap. 
This azimuthal asymmetry is at the origin of the elliptic flow in CSPM. The 
partons emitted at some point inside the cluster have to pass through the 
strong color field before appearing on the surface. The energy loss by the 
parton is proportional to the length and therefore the p t of a particle will 
depend on the direction of the emission as shown in Fig. [22] (left plot) [!158j . 
A Monte-Carlo simulation has been done taking into account this energy 
loss. The results of this simulation for different harmonics (15811159 ] are in 
reasonable agreement with the experimental data on the p t and centrality 
dependencies. In order to get an analytical expression, a simplification was 
done encoding all the azimuthal dependence in a new string density parame¬ 
ter in such a way that the inclusive azimuthal distribution is obtained by 
changing £ to £^. In this way one obtains closed formula for r; 2 (159) 116011161] 




, 2F(&)> R- 1 ' 


(130) 


The transverse momentum dependence of u 2 computed using Eq. (130) for 


Pb-Pb at yj snn = 2.76 TeV and Au-Au at = 200 GeV for 10-20% 

centrality is shown in Fig. [22] (right plot). A good agreement is also obtained 
at all centralities and rapidities \m\m- 

The ridge structure was seen first at RHIC in Au-Au and Cu-Cu and later 
at LHC in Pb-Pb collisions. This structure has been also observed in proton- 
proton and proton-nucleus collisions at high multiplicity at LHC energies. 
This strong long range rapidity correlation collimated around azimuthal an¬ 
gles 0 and 7r is one of the most interesting features of the experimental data, 
especially the structures seen in pp an pPb are a challenge for most of the 
models, in particular the hydrodynamics description. 

In simple string models azimuthal dependence was not generated, so that 
the arising correlations were flat in £>. In string percolation correlations can 
arise from the superposition of many events with different number and type 
of strings. In this way there appears long range correlations in rapidity 
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Figure 22: Left plot:Azimuthal dependence of R. R being the radius of the projected circle 
jlSSj . Right plot: Elliptic flow v 2 comparison with data and CSPM. The error bars in red 
and green are the results from Pb-Pb at yjS nn= 2.76 TeV and Au-Au at yj snn= 200 
GeV. The CSPM results are shown as dotted blue and red lines Ell¬ 


as discussed earlier. However passing to the azimuthal dependence, if the 
emission from strings is isotropic, the correlations due to their distribution 
in different events will also be isotropic. Also in the central rapidity region the 
inclusive cross section are practically independent of rapidity. This generates 
a plateau in the y — (p distribution rather than a ridge, with only a narrow 
peak at small y and 0 due to short range correlations. This conclusion 
remains true if one averages the inclusive cross sections over all events with 
the resulting loss of azimuthal angle dependence. So the ridge can only be 
obtained on an event-by-event basis. 

In string percolation the ridge arises due to fluctuations in both string 
distributions and impact parameter [lFlj . So it is important to study the 
averages relevant to the ridge formation. For a particular string configuration 
with a fixed azimuthal angle 0 O of the impact parameter, the single and 
double inclusive cross sections are 


I c (y,(p) = A c (y) + 2j2(K(y) cosra(0-0o) + C“ sinn(0-0 o )) 

n= 1 

= A c [l + 2]T(^(y)cosn(0) + <4sinn(0)) (131) 

71=1 

I c (yi,<l>i,y 2 ,(h) = U c (?/i)+2^(H^(|/ 1 )cosn(0i-0o)+C , ^sinn(0 1 -0o)) 


71=1 


X 


A c (y 2 ) + 2 ( B m(y 2 ) cos m(02 - 0o) + C^sinm(0 2 - 0o))l ■ (132) 


771=1 
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One must average this expression over string distributions and directions of 
the impact parameter. The latter reduces to integration over with weight 
1 / 2 , 71 . Doing first this integration and then averaging over string distribution 
one obtains the measured double inclusive cross section 


I(yi,y 2 ,(f>i 2 ) =< A{y 1 )A{y 2 ) > +2^ < W n (y 1 ,y 2 ) > cos (133) 

n 

where 0 12 = 0i - 0 2 and W n (y 1 ,y 2 ) = B n (y 1 )B n (y 2 ) + C n (y l )C n (y 2 ). Similar 
averaging of the single cross section eliminates all the azimuthal dependence 
I(yi) = < A(yi >. Thus the correlation function is 


C(yi,y 2 Ai2) 


< A(y 1 )A(y 2 ) > +2 E„ < W n (y u y 2 ) > cosn0i 2 
< A c > 2 


BA(y u y 2 ) 

< A(yi) >< A(y 2 ) > 


+ 2 J2 < W n(yi,y2) > COS '71012, 

n= 1 


where DA(y 1 ,y 2 ) and w n (y 1 ,y 2 ) are given by 


(134) 


DA(y u y 2 ) =< A{y-i)A{y 2 ) > - < A(y i) >< A(y 2 ) >, (135) 


w n {yi,y 2 ) 


< WJjJulte) > 
< A >2 


(136) 


The correlation contains two terms. The first is due to fluctuations in the 
total multiplicity and is independent of the angle. The ridge comes from the 
second term,which depends on the product of the initial coefficients B n and 
C n at different rapidities. 

In the central rapidity region the distributions are practically independent 
of rapidity. This creates a plateau in rapidity with the angular dependence 
of correlations given by 


C(0i 2 ) 


DA 

< A > 2 


+ 2 J2 W n COS n ^ 12 

n= 1 


(137) 


with DA = < A 2 > — < A > 2 , and w n = . If one neglects fluctuations 

of the multiplicity at a given centrality and assumes < A 2 > = < A > 2 then 
one finds w n — < 70 >, the average of the individual flow coefficients squared. 
This should be compared with the normally defined flow coefficients which 
are obtained by averaging of individual coefficients themselves. If fluctuations 
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Figure 23: Correlation function C(<j>) for (a). Au-Au at 200 GeV for 10% of the most 
central events along with experimental data (b). Pb-Pb at 2.76 TeV. (c). p-Pb collisions 
at 5.02 TeV for central collisions (d) pp collisions at 7 TeV lTi3 . 


both in A and in v n are negligible one can take < v 2 >=< v n > 2 = v 2 and 
thus find the ridge directly from the flow coefficients. 

The general scheme of the calculations of (7(0 12 ), Eq. (137), follows the 
Monte-Carlo code procedure used in the evaluations of the flow coefficients. 
In this case the color and tension of the clusters formed from n original 
strings are taken to be yfn greater than the original strings without paying 
attention to the area of the formed cluster. This simplification is done to 
save computation time to do feasible the evaluations. I 11 this way one com¬ 
putes the single and double inclusive cross-section and averaging them finds 
coefficients v n ,w n and the correlation C(0 12 ). 

Figures |23)(a) and (b) show the results [ 155] of the correlation coefficient 
C(0) for Au-Au at 200 GeV for 0-10% centrality and Pb-Pb at 2.76 TeV for 
the most central events together the experimental data [1641 [27]. A good 
agreement is observed in the regions around 0 = 0 and 7 r. It worsens at 
intermediate azimuthal angles. 

I 11 pPb collisions we assume that the maximum number of strings attached 
to the nucleons of the target to be 18 at energies in the region 5-7 TeV, in 
accordance with the results on the multiplicity on pp collisions. The number 
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of strings attached to the projectile proton will be A 1 / 3 times larger. The 
resulting w n and C(4>) are similar for energies 62.4, 200 GeV and 5.02 TeV. 
In Fig. [23[c) the correlation coefficient at 5.02 TeV for central collisions 
compared to the data [28] is shown. A very good agreement is observed. 

String percolation describes pp collisions in the same way as pA and AA 
collisions, therefore the same Monte-Carlo code procedure can be applied. 
From the multiplicity evaluations, it is known that the average number of 
formed strings at 62.4, 200 GeV and 7 TeV is 3, 4 and 9 respectively. For 
the distribution of hadronic matter in the proton a Gaussian of radius 0.8 
fin is assumed. In Fig. 23 d) we show the correlation coefficient C{<p) for 


pp collisions at 7 TeV with triple multiplicity compared to the experimental 
data [22]. We again observe a very good agreement. 




Figure 24: CSPM results on the strength of the near-side ridge for (a). Au-Au collisions 
at two RHIC energies, (b) pp collisions at 7 TeV and p-Pb collisions at 5.02 TeV 1165 1. 

Comparison with data at RHIC and LHC shows a good qualitative agree¬ 
ment in all cases. This agreement is also quantitative in cases of high multi¬ 
plicity pp and pPb collisions at 7 and 5.04 TeV respectively. The agreement 
with AA is reasonably good on the near side and on the away side but worse 
at intermediate angles. The reason for this probably has to do with simplifi¬ 
cations in the simulations done to make calculations feasible. These simplifi¬ 
cations eliminate part of the fluctuations and influence on the magnitude of 
higher harmonics. For instance, the multiplicity of a cluster was taken y/n 
instead of JnS n /S\ times the multiplicity of the original strings. Concerning 
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the pseudorapidity dependence and the strength of the near side ridge struc¬ 
ture together with the onset of it, the experimental data show an interesting 
behavior. The pp ridge structure only appears for charged multiplicity larger 
than 110 in pp collisions [29] and around 50 in pPb collisions [28] . On the 
other hand, the RHIC data on Au-Au indicate a change in the strength of 
the near side ridge structure at a given multiplicity [166]. All these facts can 
be naturally explained in string percolation [ 165] . In fact, the ridge structure 
in pp at high multiplicity was anticipated by string percolation [I67j. 

It can be shown that the strength of the ridge structure A\ defined as in 
|166] is proportional to 


dN < n(n — 1) > — < n > 2 
dy < n > 2 


(138) 


which according to Eq. (117) of the previous section becomes < N > /k 
where N = N S F (£) is the number of effective clusters and 
k =< N > /(I — exp {—£)) 3 / 2 the negative binomial parameter. Hence the 
strength of the near-side structure, A \, is proportional to 


A , = = C( 1 - exp(-0) s/2 . (139) 


The factor 1 — exp(—£) stands for an homogeneous profile. In the case of 
more realistic profiles as Gaussian or Wood-Saxon we should use Eq. (15), 
thus 

Ai = C- --- 7372 , ( 14 °) 


(l + aexp(-(f-&))/& 


where C is a normalization constant which is different for each type of colli¬ 
sion. 

Comparison of the experimental data for the dependence of Au-Au at 
62.4 and 200 GeV on the multiplicity [ 166J and pPb and pp at 5.04 and 7 
TeV respectively [28] is shown in Figs. [24|a) and (b). In these figures the 
parameters a and £ c are kept fixed, with values 1.5 and 1.5 respectively (with 
these values one obtains the correct dependence of the fraction of the area 
covered by strings on the string density). The values of b are 0.75 for Au-Au 
collisions and 0.35 for both pPb and pp collisions. We would expect b ~ 1/R, 
with R the transverse dimension of the collision area, which is not far from 
the values found. Also a good agreement is found [165] for the behavior 


65 

















of the width in pseudorapidity of the ridge structure as a function of the 
multiplicity in Au-Au collisions at 62.4 and 200 GeV. 

We conclude that string percolation is able to describe the ridge struc¬ 
tures seen in pp high multiplicity, pPb and AA collisions. In this framework 
the ridge is obtained from the superposition of many events with different 
numbers and types of clusters of strings. In order to obtain anisotropy it 
is crucial to introduce quenching of the emitted partons in the strong color 
held inside the clusters of strings. In this way, it can be seen that the ridge 
structure in produced by an interplay of initial stage effects (formation of 
different clusters) and final state effects described by quenching. 

3-4- Bose-Einstein correlations (BEC) 

Studies of correlations between two identical particles have been very 
fruitful to determine the extension of the source of multiparticle production 
in collisions of two hadrons. This has been used extensively in heavy ion 
experiments at RHIC at LHC |1681 11691 11701 11711 11721 11731 1174] . These 
correlations provide information not only on the geometrical structure but 
also on the degree of coherence of the emitted particles. The correlation 
strength is characterized by a parameter A, which can also be interpreted as 
a measure of chaoticity or the degree of coherence of the collision [1751 11761 
1177] , In this interpretation A = 1 signals totally chaotic emission whereas A 
= 0 means radiation in a coherent way. However, this interpretation should 
be taken with caution. In fact, in e + — e~ annihilation, A is close to one 
at energies where there is no production of three or more jets and in the 
more complex pp collisions A decreases with increasing multiplicity. These 
observations would apparently suggests a systematic increase of the coherence 
from e + — e~ to pp collisions, which does not seem reasonable. 

There are many difficulties in extracting the value of A from measure¬ 
ments (disentanglement from the Coulomb repulsion, corrections from long 
lived resonances, extrapolation to low p t , corrections due to experimental 
transverse resolution). However, the existing measurements at ISR, SPS, 
RHIC and LHC give valuable information which should be understood. The 
experimental data on A have been obtained in different kinematic situations 
and also assuming different extrapolations, corrections and parametrizations, 
which makes rather difficult their comparison with theoretical models or even 
with different data. However some trends can be distinguished. First, for a 
not very large number of collisions, the data at SPS with p and Oxygen pro¬ 
jectiles show a decrease of A with multiplicity [1781117911180] , As the number 


66 























of collisions increases the behavior of A changes. It no longer decreases with 
multiplicity and even may increase. This is true at LHC where A reaches the 
values of 0.7-0.74. At SPS energies, the A values found are larger at forward 
than at central rapidity (Notice that the particle multiplicity is larger at cen¬ 
tral than in forward rapidity). This behavior can be naturally understood in 
the framework of percolation of strings [IBTl, [182] . 

The strings of the Lund type correspond to totally chaotic sources [1831 
1184] . A = 1, and usually it is assumed that there is no BEC from a pair of 
particles emitted from different strings (a discussion on this assumption can 
be seen in m\)- Under these assumption one finds [ lEtJ j 

A = — (141) 

T^tot 


where n s is the mean number of identical particle pairs produced by frag¬ 
mentation of the same string in a given rapidity and transverse momentum 
range and n tot is the total number of identical pairs produced in the same 
kinematic range. The number of pairs of identical particles produced by each 
cluster is 


n s 


1 2 v—' 

n/U < L “AT" 

Z n =1 


> 


(142) 


and the total number of pair of identical particles produced is 

2 


ntot = 2^1 < 


'^ a n y/nS n 


Kn =1 


Si 


> 


(143) 


where a n is the number of clusters with n strings, number of particles pro¬ 
duced by a single string in the considered range. A Monte-Carlo simulation 
taking into account energy conservation was done at several energies and dif¬ 
ferent projectiles and targets [181] . It was found that there is an approximate 
scaling in the string density, i.e. the chaoticity A only depends on it. In Fig. 


25 values of A are shown as a function of the string density £. The details of 


the Monte-Carlo calculations together the SPS experimental results and the 
predictions for RHIC and LHC energies are shown in Ref. [ 181 ] . We observe 
a good agreement with the data in particular with the change of the behavior 
with multiplicity. Notice that the value of A = 0.7 predicted at LHC (last 
point in Fig. (25)) was confirmed by the LHC data [ 172 . 173, 174 ] . The three 
body BEC have been also studied in the framework of percolation [182], pre¬ 
dicting the strength of the genuine three particle Bose-Einstein correlations 
and obtaining a general agreement with experimental data [1741118711188] . 
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Figure 25: Dependence of A on £ for different nucleus-nucleus collisions with percolating 
strings 11811 . 

3.5. .//'I' Suppression 

Melting of the J/\k due to the Debye screening in a deconhned medium 
was one of the first proposed signatures of the existence of quark gluon plasma 
|i|. Very early, at SPS energies, it was found that J/'k was suppressed 
in S-U and O-U collisions [ 190 . 1 191 ] . However this suppression was naturally 
explained by absorption of the pre-resonant cc state in the nucleus [ 1921 193]. 
From that time, many efforts have been devoted to disentangle the suppres¬ 
sion due to nuclear effects from the actual melting in a deconhned medium. 
The RHIC and the LHC experimental data pointed out that, indeed, in ad¬ 
dition to the nuclear effects there is a sequential suppression of J/T [194] 
combined with an enhancement due to a recombination mechanism ra 
and to the existence of strong color fields above a certain matter density. On 
the other hand, ALICE collaboration [ 196 ] have shown that J/T production 
in pp high multiplicity events, normalized to the minimum bias, as a func¬ 
tion of the multiplicity in the central rapidity region also normalized to the 
minimum bias increases linearly or even more rapidly at high multiplicity. 

In string percolation it is expected that above the critical percolation 
threshold the J/T or at least the T's are suppressed. However, the string 
density reached at 7 TeV is much smaller than the percolation critical value 
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Figure 26: Our results for pp collisions in the central (|y| < 0.9, dashed line) and forward 
(2.5 < y <4, dotted line) rapidity range, together with the experimental data for the 
central (circles) and forward (squares) rapidity regions from the ALICE collaboration 
1196] . The linear behavior(solid line) and our prediction for p-Pb collisions (solid line) at 
7 TeV are also plotted [197] . 
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even for high multiplicity events reached in the experiment. The behavior 
found by the ALICE collaboration can be naturally explained in the frame¬ 
work of string percolation |197j . In fact, assuming that, as in any hard pro¬ 
cess, the number of produced J/T is proportional to the number of collisions 
and thus is proportional to the number of strings N s , we can write 


nj/m _ N s 
< nj/m > ~ < N s >' 

On the other hand we have 

dN/dy N a F(£) 

< dN/dy > <N S > F(< £ >)' 

From above equations we have 

dN/dy = / nj/m \ 1/2 ( 1 ~ ex P(~ <Zj> \ 
< dN/dy > y< nj/^ > ) I 1 — exp(— < £ >) 1 


(144) 


(145) 


(146) 


At low multiplicities where the number of strings N s is small, the above 
equation gives rise to the linear dependence 

rij/xy dN/ dy 

< dN/dy > 


(147) 


< nj/xs, > 

On the othe r hand, at high multiplicities, the bracket on the right hand side 
of Eq. (146) can be approximated by < £ > -1 / 2 . One obtains in this case 


nj /\Jr 


=<£> 


dN/dy 
< dN/dy > 


(148) 


< nj/* > 

Thus the linear dependence obtained previously changes to quadratic when 
the high multiplicity events are in play. At much higher multiplicity, about 
7-8 times minimum bias we predict a rapid suppression of this ratio. At 14 
TeV the string density in pp increases by about 20% and thus at multiplicities 
of 5-7 times greater than at minimum bias there will be a suppression of J/T 
free of nuclear effects. 

In Fig. (26) we show the results from Eq. (146) (dashed line) in the 
central rapidity region together the experimental data |196j (circles). Also the 
results for the forward region (dotted line) together with the data (squares) 
are shown. A very good agreement is observed. Note that in the forward 
direction there are less number of strings, producing less multiplicity and 
thus the departure of the linear behavior would be at higher multiplicities. 
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4. Thermodynamic and transport properties 

A possible phase transition of strongly interacting matter from hadron 
to quark-glnon plasma state has received considerable interest in the past. 
What conditions are necessary in the pre-equilibrium stage to achieve de¬ 
confinement and perhaps subsequent Quark Gluon Plasma(QGP) formation 
? 

As discussed earlier in Section 2.2 the multiplicity /j and the mean trans¬ 
verse momentum squared (pf) of the particles produced by a cluster of n 
strings |198| 

Pn = (Pt) = i ( 149 ) 

where /ii and (pl)i are the mean multiplicity and (pi) of particles produced 
from a single string with a transverse area S\ = nr q. In the thermodynamic 
limit, one obtains an analytic expression [T7J 1198] 


<^> = 


1 - e-f F(if 


m 



(150) 


Here £ = is the percolation density parameter with Ng and Sn being 
the total number of strings and interaction area respectively. F(£) is the 
color suppression factor and is related to a n and (pf) n through the relation 
Pn = F(£)/x 0, ( pl)n = (p# 2 )i/F(£). 

In this section thermodynamic and transport properties are discussed in 
the framework of clustering of color sources at RHIC and LHC energies. The 
thermodynamics of clustering can be addressed by extracting the temper¬ 
ature from the transverse momentum spectra of charged hadrons. In the 
model a local temperature can be introduced. In fact the tension of a cluster 
of strings fluctuates around its mean value because the chromoelectric held 
is not universal. Assuming a Gaussian form for these fluctuations a thermal 
spectrum is obtained with a temperature inversely proportional to ^2F(£). 
It is remarkable that for £ = 1.2 — 1.5 the corresponding temperature T c is 
not far from the critical temperature [ 199] . We extract the color suppression 
factor F(£) using experimental data at RHIC energies. F(£) and £ permit 
the determination of the initial temperature T and the Bjorken initial en¬ 
ergy density e. Energy density e was found to be a linear function of £. As a 
result one can derive the equation of state and the transport coefficient shear 
viscosity over entropy density ratio as a function of the temperature. A good 
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agreement between the model and lattice QCD results for energy density is 
obtained for 1 < T/T c < 2 [200] . Relativistic kinetic theory expression for 
shear viscosity has been used to obtain shear viscosity over entropy ratio rj/s 
at RHIC and LHC energies [201J. 77 /s has a minimum at the temperature 
close to the critical one, rising slowly as the temperature increases. The 
data based values are compared with theoretical values for both a weakly 
coupled QGP (wQGP) and a strongly coupled QGP (sQGP). The dimen¬ 
sionless quantities, the trace anomaly A = (e — 3 p)/T 4 and the inverse of 
77 /s, have identical magnitudes for 0.5 < T/T c < 6 which describes the tran¬ 
sition from wQGP to sQGP and then back to wQGP. This ’Ansatz’ provides 
an independent value for pressure at a given e | 2 U 2 j . The sound velocity 
Cf obtained using this pressure is in excellent agreement with the LQCD 
simulations [ 200 ]. 

4-1. Experimeiital determination of the color suppression factor F(£) 

To evaluate the initial value of £ from data, a parameterization of p-p 
events at 200 GeV is used to compute the p t distribution |203j 



Figure 27: Charged particle pt spectrum from pp collisions (blue color) and Au-Au 
collisions at 200 GeV for 0-10% centrality |109| . 


dNJdpl = a/(p 0 +PtY 


(151) 


where a, p 0 , and a are parameters used to fit the data. Fig. 27 shows the 


normalized pt spectrum for pp collisions at a/s = 200 GeV. A fit to the 
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data gives p 0 = 1.982 and a = 12.877. This parameterization is used for 
nucleus-nucleus collisions to take into account the interactions of the strings 


[293] 


/ (nSi/S n ) Au ^ 

V <nSi/S n > 


1/4 


(152) 


where S n corresponds to the area occupied by the n overlapping strings. The 
thermodynamic limit, i.e. letting n and S n —>■ oo while keeping £ fixed, is 
used to evaluate 


nS\ 

~S 


i/f 2 K) 


(153) 


dN c /dpi 


a 

Wnw/no+Pif 


(154) 


In pp collisions ( nSi/S n ) pp ~ 1 due to the low string overlap probability. 
The STAR analysis of charged hadrons for Au-Au collisions at yfs/vjv = 200 
and 62 GeV for £ is shown in Fig. [28] as a function of number of participants 


[203]. It is observed that for all centralities, except for the most peripheral 
one, yj snn = 200 GeV data lie above the critical percolation threshold £ c = 
1.2. The value of £ increases with increasing collision centrality, an expected 
indication of additional string overlap in more central collisions. For y/s^N 
= 62.4 GeV, almost all centralities, except for the three most peripheral one 
lie above £ c . 

Figure 29 shows a plot of F(£) as a function of charged particle multi¬ 
plicity per unit transverse area for Au-Au collisions at yj snn = 200 

GeV for various centralities for the STAR data [10911199] . The error on F(£) 
is ~ 3%. R(£) decreases in going from peripheral to central collisions. The 
£ value is obtained using Eq. (150), which increases with the increase in 
centrality. The fit to the Au-Au points has the functional form 


F(£) = exp[—0.165 — 0.094—^—^ / S' at] - 


(155) 


Recently, the ALICE experiment at LHC published the charged-particle mul¬ 
tiplicity density data as a function of centrality in Pb-Pb collisions at y / Sjv r /v 
= 2.76 TeV m The STAR results for Au-Au collisions at = 200 

GeV have been used to estimate F(£) values for Pb-Pb collisions at different 


centralities using the fit function Eq. (155) for Au-Au. 
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Figure 28: Percolation density parameter £ as a function of number of participants for 
Au-Au collisions at yj snn = 62 and 200 GeV [203j . 



the size of the symbol. The line is fit to the STAR data. The solid blue 


. The solid red 
is smaller than 
squares are for 


Pb-Pb at 2.76 TeV [201]. 
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4-2. Energy density 

The QGP according to CSPM is born in local thermal equilibrium because 
the temperature is determined at the string level. After the initial tempera¬ 
ture T > T c the CSPM perfect fluid may expand according to Bjorken boost 
invariant ID hydrodynamics [7j 



dNc 

dy 


(m t ) 


Sn Tpro 


(156) 


where £ is the energy density, S n nuclear overlap area, and T pro the production 
time for a boson (gluon) [204] . 

The dynamics of massless particle production has been studied in QE2 
quantum electrodynamics. QE2 can be scaled from electrodynamics to quan¬ 
tum chromodynamics using the ratio of the coupling constants [204] . The 
T pro for a boson (gluon) is [205] 


_ 2.405ft 

{m t ) 


(157) 


Above the critical temperature only massless particles are present in 
CSPM. To evaluate £ we use the charged pion multiplicity dN c /dy at midra¬ 
pidity and S n values from STAR for 0-10% central Au-Au collisions with 
\J$nn =200 GeV [199]. The factor 3/2 in Eq. (156) accounts for the neutral 

pions. The average transverse mass (m t ) is given by (m t ) = \J ( p t ) 2 + Wq, 
where ( p t ) is the transverse momentum of pion and m 0 being the mass of 
pion . 

From the measured value of £ and e, as shown in Fig. 

£ is proportional to £ for the range 1.2 < £ < 2i 


30 


it is found that 
£ = 0.788 £ (GeV//m 3 ) 


[ 201] . The extrapolated value of £ for central Pb-Pb collision at 2.76 TeV is 
8.32 GeV/ fm 3 as shown in Fig. 


30 


4-3. Temperature measurement and thermalization 

The connection between £ and the temperature T(£) involves the Schwinger 
mechanism (SM) for particle production. The Schwinger distribution for 
massless particles is expressed in terms of p\ 205 . 1204] 


dn/dpl ~ exp(— irp/fx 2 ) 


(158) 
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Figure 30: Energy density e as a function of the percolation density parameter £. The 
extrapolated value for LHC energy is shown as blue square BO 11 . 


where the average value of the string tension is ( x 2 ). The tension of the 
macroscopic cluster fluctuates around its mean value because the chromo¬ 
electric field is not constant. The origin of the string fluctuation is related 
to the stochastic picture of the QCD vacuum. Since the average value of 
the color field strength must vanish, it can not be constant but changes 
randomly from point to point [rrcnrcnfi] . Such fluctuations lead to a Gaussian 
distribution of the string tension 


dn 

dpt 



~2 2 
)exp(—7M 


2 < x 2 > 


x- 


which gives rise to thermal distribution [103] 


(159) 


dn . I 27r 

«~ exp( - p *VF) ) ' 


(160) 


with ( x 2 ) = 7r(pt)i/F(£). The temperature is expressed as [1041 1199] 


no 


Mi 

\ 2F(f)' 


(161) 
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We will adopt the point of view that the experimentally determined chem¬ 
ical freeze-out temperature is a good measure of the phase transition tem¬ 
perature, T c [207] . The single string average transverse momentum (p 2 )i is 
calculated at £ c = 1.2 with the universal chemical freeze-out temperature of 
167.7 ± 2.6 MeV [208]. This gives sj(p 2 )i = 207.2 ± 3.3 MeV which is close 
to ~ 200 MeV used previously in the calculation of percolation transition 
temperature ES- 

Recently, it has been suggested that fast thermalization in heavy ion 
collisions can occur through the existence of an event horizon caused by a 
rapid deceleration of the colliding nuclei [ 209] . The thermalization in this 
case is due the Hawking-Unruh effect |210l 1211] , 

It is well known that the black holes evaporates by quantum pair produc¬ 
tion and behave as if they have an effective temperature of 


T h = 


1 

8ttGM’ 


(162) 


where 1/4GM is the acceleration of gravity at the surface of a black hole of 
mass M. The rate of pair production in the gravitational background of the 
black hole can be evaluated by considering the tunneling through the event 
horizon. Unruh showed that a similar effect arises in a uniformly accelerated 
frame, where an observer detects the thermal radiation with the temperature 
T =a/2, where a is the acceleration. Similarly, in hadronic interactions the 
probability to produce states of masses M due to the chromoelectric held E 
and color charge is given by the Schwinger mechanism 

— 7T M 2 

Wm ~ exp(-——) ~ exp (-M/T), (163) 

gE 


which is similar to the Boltzmann weight in a heat bath with an effective 
temperature 


T=— ,a= 2 ^. (164) 

2t r ’ Af V ; 

In CSPM the strong color held inside the large cluster produces de¬ 
acceleration of the primary qq pair which can be seen as a thermal tem¬ 
perature by means of Hawking-Unruh effect. This implies that the radia¬ 
tion temperature is determined by the transverse extension of the color hux 
tube/cluster in terms of the string tension [ 2121 213 ]. 


T = 



(165) 
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The string percolation density parameter £ which characterizes the per¬ 
colation clusters measures the initial temperature of the system. Since this 
cluster covers most of the interaction area, this temperature becomes a global 
temperature determined by the string density. In this way at £ c = 1.2 the 
connectivity percolation transition at T(£ c ) models the thermal deconfine¬ 
ment transition. The temperature obtained using Eq. (161) was ~ 193.6 



Figure 31: Average < pt > of pions as a function of temperature obtained using Eq. (161). 


MeV for Au-Au collisions at yj snn = 200 GeV in reasonable agreement with 
T t = 221 ±19 stat ± 19 sys from the enhanced direct photon experiment mea¬ 
sured by the PHENIX Collaboration [IS]. For Pb-Pb collisions at y / sjv r /v 
= 2.76 TeV the temperature is ~ 262.2 MeV for 0-5% centrality, which is 
expected to be ~ 35 % higher than the temperature from Au-Au collisions 
[ 199] . A recent summary of the results from Pb-Pb collisions at the LHC 
has mentioned that the initial temperature increases at least by 30 % as 
compared to the top RHIC energy m- The direct photon measurements 
from ALICE gives the temperature of T % = 304 ±51 |215j . The agreement 
with the measured temperaure shows that the temperature obtained using 
Eq. (161) can be termed as the initial temperature of the percolation cluster. 
For an ultrarelativistic ideal gas ( p t ) oc T. Figure [3l] shows a plot of ( p t ) as a 
function of the temperature. When the transverse momentum is exponetially 
distributed with inverse slope T in a given event, ( p t ) = 2 T [216] . Using 
this relation one can predict the ( p t ) for pions in Pb-Pb collisions at y/SNN = 
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Table 1: The percolation density parameter temperature T, rp?’o, and y/s for meson 
gas [2H], hadron to QGP transition in Au-Au collisions at 200 GeV and Pb-Pb at 2.76 
TeV (estimated). Au-Au is for 0-10% and Pb-Pb is for 0-5% central events. 


System 

£ 

T (MeV) 

v/s 

1 

O 

s- 

G- 

Meson Gas 

0.22 

150.0 

0.76 


Hadron to QGP 

1.2 

167.7 ±2.6 

0.240T0.012 


Au-Au 

2.88T0.09 

193.6 ±3.0 

0.204±0.020 

1.13 

Pb-Pb 

10.56il.05 

262.2±13.0 

0.260T0.026 

0.89 


2.76 TeV. For 0-5 % the (p t ) = 0.53 GeV/c in agreement with the measured 
value of 0.517 ± 0.019 GeV/c |217] . Table I gives the CSPM values £, T, 
and e at T/T c = 0.88, 1, 1.16 and 1.57. 
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Figure 32: e/T 4 versus T/T c from CSPM (red circles) and Lattice QCD (blue dash line) 
for 2+1 flavor and p4 action [TW| . 


4-4- Shear viscosity 

The observation of the large elliptic flow at RHIC in non-central heavy 
ion collisions suggest that the matter created is a nearly perfect fluid with 
very low shear viscosity nansiininBi. Recently, attention has been focused 
on the shear viscosity to entropy density ratio rj/s as a measure of the flu¬ 
idity (2191122011221L 222 j. The observed temperature averaged rj/s based on 
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viscous hydrodynamics analyses of RHIC data , are suggestive of a strongly 
coupled plasma |22311224] . The effect of bulk viscosity is expected to be neg¬ 
ligible. It has been conjectured, based on infinitely coupled super-symmetric 
Yang-Mills (SYM) gauge theory using the correspondence between Anti de- 
Sitter(AdS) space and conformal field theory (CFT), that the lower bound 
for r)/s is 1 / 47 T and is the universal minimal r\js even for QCD [19]. However, 
there are theories in which this lower bound can be violated [225]. 



Figure 33: A m fp in fm ( blue line), temperature in GeV scaled by a factor of 10(green 
line) and A m / p x scaled T in GeV fm (red line) as a function of £. The minimum in 77 /s is 
due to the combination of A m f p and T. Meson gas [218j . RHIC and LHC points are also 
shown as solid black square 1201] . 


The relativistic kinetic theory relation for the shear viscosity over entropy 
density ratio, rj/s is given by [22311226] 


] ^ T\ m f p 

s ~ 5 


(166) 


where T is the temperature and A m fv i s the mean free path. A m f P ~ (TA 77 
where n is the number density of an ideal gas of quarks and gluons and a tr 
the transport cross section. 

I 11 CSPM the number density is given by the effective number of sources 
per unit volume [133, 227] 


n = 


N„. 


SnL 


(167) 
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L is the longitudinal extension of the string ~ 1 fm [104] . The area occupied 
by the strings is related to the percolation density parameter £ through the 
relation (1 — e _ ^)S/v- Thus the effective no. of sources is given by the total 
area occupied by the strings divided by the area of an effective string SiF(£). 




(1 - e~*)S N 
SiF(t) 


(168) 


In general N. 


sources 


is smaller than the number of single strings. N s , 


equals the number of strings N s in the limit of £ = 0. The number density 
of sources from Eqs. (167) and (168) becomes 


(1 ~ e~ € ) 

SiF(0L ' 


(169) 


The transport cross section cr tr is the transverse area of the effective string 
S\F(£). Thus <jfr is directly proportional to A, which is in agreement with 
the estimated dependence of a tr on the temperature [ 2261 228] . The mean 
free path is given by 

Xmfp = (i_e-*)‘ 

For large value of £ the A m f p reaches a constant value. r)/s is obtained from 
£ and the temperature 

77 _ TL 
s ~ 5(1 -e-f) 

Below £ c , as the temperature increases, the string density increases and the 
area is filled rapidly and X m f P and rj/s decrease sharply. Above £ c , more than 
2/3 of the area are already covered by strings, and therefore the area is not 
filling as fast and the relatively small decrease of A m j p is compensated by the 
rising of temperature, resulting in a smooth increase of 77 /s. The behavior 
of 77 /s is dominated by the fractional area covered by strings. This is not 
surprising because 77 /s is the ability to transport momenta at large distances 
and that has to do with the density of voids in the matter. Figure [33] shows a 
plot of A m f p , T and A m f p x T as a function of £. Thus the product T(£) x \ m f p 
will have a minimum in 77 /s. It has been shown that 77 /s has a minimum 
at the critical point for various substances for example helium, nitrogen and 
water [ 224] . Thus the measurement of 77 /s as a function of temperature 
can locate the critical end point/crossover region in the QCD phase diagram 
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Figure 34: r]/s as a function of T /T c . Au-Au at 200 GeV for 0-10% centrality is shown 
as solid black square. wQGP and sQGP values are shown as dotted blue and green lines 
respectively 13531 . The red dotted line represents the extrapolation to higher temperatures 
from the CSPM [201 . The hadron gas value for 77 /s ~ 0.7 is shown as solid black circle 
at T /T c ~0.88 |218: . The red solid line is Ads/CFT limit [T9l . 



T ~ 175-185 MeV. Figure 34 shows a plot of 77 /s as a function of T/T c . 
The estimated value of rj/s for Pb-Pb is also shown in Fig. 34 at T /T c = 

The 


1.57. The lower bound shown in Fig. 34 is given by AdS/CFT 
results from STAR and ALICE show that the 77 /s value is 2.5 and 3.3 times 
the KSS bound [19]. Extrapolating the 77 /s CSPM values it is clear that it 
will approach the weak coupling limit near T/T c ~ 5.8 as the temperature is 
increased. The CSPM 77 /s value for the hadron gas is in agreement with the 
calculated value using measured elastic cross sections for a gas of pions and 
kaons [218] . 


4 .5. Jet quenching parameter 

One of the outstanding observation made in central collisions of heavy 
ions is that the final state hadrons are strongly suppressed at large trans¬ 
verse momenta This phenomena is referred as jet quenching, is 

understood to be caused by gluon radiation induced by multiple collisions of 
the leading parton while traversing the medium [231] , In QCD, the energy 
loss is expressed by the transport coefficient q, which describes the rate at 
which medium transfers transverse momentum of the fast parton. An impor- 
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tant question is whether different transport coefficients are related with each 
other. It was shown that two transport coefficients shear viscosity 77 and jet 
quenching parameter q are related [231, l232j. 


n 3 T 3 
_ _ 

s 2 q 


(172) 


The relation associates a small ratio of rj/s to a large value of q. A large 
amount of theoretical work has been done to extract the jet transport pa¬ 
rameter from jet quenching at RHIC and LHC energies [229, 1221 - [2221 222. 
The latest study by the JET Collaboration has extracted or calculated q 
from five different approaches to the parton energy loss in a dense medium. 
The evolution of bulk medium in the study was given by 2+1D or 3+1D hy¬ 
drodynamic models with the initial temperatures of Tr#/™ = 346 - 373 MeV 
and T^q 0 = 447 - 486 MeV for most central Au-Au collisions at i/sjviv = 
200 GeV and Pb-Pb collisions at i/sjviv = 2.76 TeV respectively. The vari- 
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Figure 35: Scaled jet quenching parameter p as a function of the temperature. The 
values shown in solid blue squares (Hydro(RHIC) and Hydro(LHC) are given by Eq. (173). 
The CSPM values are shown in green solid squares CSPM(RHIC) and CSPM(LHC) for 
temperatures ~193 and ~262 MeV at RHIC and LHC energies respectively [20l j. 


ation of q values between different models can be considered as theoretical 
uncertainties. One therefore can extract its range of values at RHIC and 


83 


















LHC [12251 233 ]. The scaled q at the highest temperatures reached in the 
most central Au-Au collisions at RHIC and Pb-Pb collisions at LHC are 


( f ^ r4.5±1.3 at RHIC 
rj -\3 r ''-' I3.7d=1.4 at LHC * 


(173) 


The corresponding absolute values for q(GeV 2 / fm) for a 10 GeV quark jet 
are, 

rl.2±0.3 T=370MeV 


q ~ {l 


9±0.7 T=470MeV> 


(174) 

at an initial time r 0 = 0.6 fm/c. The temperature dependence of scaled jet 


transport parameter is shown in Fig. 35 The CSPM values for using 


the relation given by Eq. (172), are shown in Fig. 35 (solid red circles). The 
CSPM values for RHIC and LHC energies are shown as solid green squares 
while the theoretical values are shown as blue squares. It is observed that 
CSPM values are in agreement with the JET Collaboration results. 

The phase transition and q have been studied in the framework of dy¬ 
namical holographic QCD model |234j . It is found that both q/T 3 and trace 
anomaly A = (e — 3 p) show a peak around the critical temperature T c . This 
indicates that q can characterize the phase transition [234]. 

In the next section we show that in CSPM the inverse of 77 /s is related 
to the trace anomaly. 

4-6. Trace anomaly 

The trace anomaly (A) is the expectation value of the trace of the energy- 
momentum tensor, ( 0 /) — (e — 3 p), which measures the deviation from 
conformal behavior and thus identifies the interaction still present in the 
medium [235]. We find that the reciprocal of rj/s is in quantitative agreement 


with (e — 3 p)/T A over a wide range of temperatures [ 2021 236]. This result 
is shown in Fig. 36 The minimum in 77 /s ~ 0.20 determines the peak of the 
interaction measure ~ 5 in agreement with the recent HotQCD values [ 200] , 


This happens at the critical temperature of T c ~ 175 MeV. Figure 36 also 
shows the results from Wuppertal Collaboration [237]. The maximum in A 
corresponds to the minimum in 77 /s. Both A and 77 /s describe the transition 
from a strongly coupled QGP to a weakly coupled QGP. 


4-7. Equation of State: Sound velocity 

We use CSPM coupled to a ID Bjorken expansion. The input parameters 
are the initial temperature T, the initial energy density e, and the trace 
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Figure 36: The trace anomaly A = (e — 3p)/T 4 vs temperature [200 y Red solid squares are 
from HotQCD Collaboration [200] . Black stars are from Wuppertal Collaboration [237] . 


anomaly A are determined by data. The Bjorken ID expansion gives the 
sound velocity 


1 dT 

T~Tr 

dT 

dr 

de 

dr 


~C 2 Jr 

(175) 

dT de 

(176) 

de dr 

~Ts/t , 

(177) 


where £ is the energy density, s the entropy density, r the proper time, and 
C s the sound velocity. Since s = (e + P)/T and P = (e — AT 4 )/3 one gets 


& 1 > = < 178 ) 

From above equations Cl can be expressed in terms of £ 

C = (-0.33) (ACh _ +0.0191 (A/3) ( (1 ^ )2 - (179) 


Since there is no direct way to obtain pressure in CSPM, we have made the 
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Figure 37: The speed of sound from CSPM versus T/T ^ SPM (red circles) and Lattice 
QCD-p4 speed of sound versus T/T^ CD (blue dash line) [238] . 
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assumption that A = (e — 3 P) ~ 1 /( 77 /s) ■ Fig. 37 shows a plot of as 
a function of T/T c . It is observed that the CSPM results are in very good 
agreement with the lattice calculations [238J. This suggest that the A can 
be approximated to l/{rj/s). 


5. Summary 

The string percolation gives a successfully description of most of the ex¬ 
perimental data in the soft region, including rapidity distributions, proba¬ 
bility distributions of multiplicities and transverse momentum, strength of 
the B-E correlations, elliptic flow and ridge structures. The clustering of 
color sources has a clear physical basis, although is not derived directly from 
QCD, but has its most fundamental features. The non-Abelian character is 
reflected in the coherence sum of the color fields which as a consequence gives 
rise to an enhancement of the transverse momentum and a relative multiplic¬ 
ity suppression. The conhnement is reflected in the small transverse size of 
the strings as well as in the length of the transverse correlations. The scaling 
observed in the transverse momentum distribution is a consequence of the 
invariance of the probability distributions under the transformations of the 
size of the clusters of strings. 








The interpretation of data using string CSPM describes the equation of 
state close to the critical temperature and also obtains a reasonable depen¬ 
dence on the temperature of the shear viscosity over entropy density ratio. 
Thermodynamical results are also in agreement with the lattice QCD calcu¬ 
lations. The percolation critical transition is known to represent a continuous 
phase transition. In central Au-Au collisions at y / Iw/v = 200 GeV the QCD 
to hadron phase transition for baryon density hb ~ 0 is believed to be a 
cross-over transition which does not have a latent heat. The CSPM EOS 
correctly describes the QCD to hadron cross-over transition and provides an 
answer to the question of the origin of universal temperature observed in 
A-A, p-p and e + e _ collisions. The percolation framework provides us with 
a microscopic picture which predicts the early thermalization required for 
hydrodynamical calculations. 

The minimum in 77 /s can be studied as a function of the beam energy 
at RHIC that could locate the critical point/crossover in the QCD phase 
diagram seen in substances like helium, nitrogen and water. The accurate 
determination of 77 /s is also important for the evaluation of another transport 
coefficient, the jet quenching parameter q. 

There are many similarities between string percolation and the Glasma 
of the color glass condensate, in formal aspects and also in phenomenological 
applications. In this way string percolation can be regarded as a complemen¬ 
tary view of CGC of the initial stage of a high hadronic collision. 
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